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1.0  INTRODUCTION 

With  the  availability  of  advanced  computer  codes,  such  as  the  Dynamics  of 
Rolling  Element  Bearings  (DREB)  computer  program  [1-5],  the  interest  in 
real  time  performance  simulation  of  rolling  bearings  has  been  growing 
rather  rapidly  over  the  past  few  years.  The  most  generalized  dynamic 
model  underlining  DREB  has  been  employed  for  the  purpose  of  both  design  and 
performance  diagnosis  in  advanced  rolling  bearing  systems.  However,  in 
spite  of  the  rapidly  advancing  computer  technology,  the  computational  effort 
required  in  using  DREB  has  proven  to  be  a  practical  limitation.  This  is 
related  directly  to  the  nature  of  the  model  in  the  sense  that  the  maximum 
step  size  for  the  numerical  integration  of  a  system  of  differential  equations 
of  motion  is  determined  by  the  highest  frequency  present  in  the  system,  which 
is  often  several  orders  of  magnitude  higher  than  the  shaft  rotational  speed. 
Thus  in  order  to  simulate  effects  of  the  order  of  shaft  speed,  integration 
over  a  very  large  time  domain  is  required  and  since  the  maximum  step  size 
is  limited  by  the  high  frequencies,  the  computing  costs,  specially  for  very 
low  speed  applications,  have  proven  to  be  prohibitive.  The  primary  objective 
of  the  present  investigation  is  to  examine  the  various  time  scales  in  the 
general  behavior  of  a  ball  bearing  and  accordingly  modify  DREB  to  suppress 
the  very  high  frequencies  and  thereby  develop  a  "rapid  version"  of  DREB, 
which  has  been  called  RAPIDREB. 

For  a  typical  ball  bearing  the  high  frequency  present  in  the  system  is 
defined  by  the  ball/race  Hertzian  contact.  If  the  ball  is  free  to  translate, 
vibratory  motion  corresponding  to  this  natural  frequency  has  been  observed. 
Also,  as  will  be  discussed  in  the  next  section,  a  relatively  high  frequency 
is  also  associated  with  the  kinematics  of  ball  motion.  In  order  to  suppress 
these  high  frequencies  the  ball  mass  center  may  be  constrained  in  a  suit¬ 
able  manner.  Once  this  is  done  the  step  size  can  be  made  relatively  large 
and  also  predictor-corrector  type  algorithms  can  be  employed  to  consider¬ 
ably  reduce  the  computing  costs.  However,  if  a  cage  is  present 
in  the  bearing,  the  discontinuities  associated  with  the  ball/cage  and  race/ 
cage  collisions  can  produce  stability  problems  for  the  predictor-corrector 
schemes  and  therefore  such  algorithms  can  only  be  effectively  used  in  the 
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case  of  a  cageless  bearing,  and  considerable  care  will  be  required  in  using 
such  algorithms  when  there  are  repeated  ball/cage  or  race/cage  collisions. 

In  general  it  will  be  necessary  to  re-start  the  predictor-corrector  process 
at  the  time  of  the  start  and  end  of  such  collisions  if  the  collision  time 
is  large.  For  short  collisions  an  explicit  formula  can  be  used  during  the 
entire  collision.  The  ultimate  integrating  scheme  may  therefore  be  of  a 
hybrid  nature  employing  both  explicit  and  implicit  algorithm.  Thus,  the 
incorporation  of  both  types  of  formulae  has  been  one  of  the  objectives  of  the 
present  work. 

The  next  section  of  this  report  reviews  the  various  time  scales  associated 
with  generalized  motion  of  the  rolling  elements  in  a  ball  bearing.  Section 
3  contains  the  essential  modifications  to  the  existing  program  DREB  in  order 
to  develop  the  RAPIDREB  version.  Some  dynamic  simulations  for  both  the  very 
low  speed  DMA*  bearings  and  high  speed  engine  bearings  are  then  presented  as 
typical  test  cases  demonstrating  the  capabilities  of  RAPIDREB.  Descriptions 
of  program  inputs  are  presented  in  Appendix  A,  and  typical  outputs  of  the 
program  are  shown  in  Appendices  B  and  C. 


*  Despun  Mechanical  Assembly 


2.0  GENERAL  TIME  SCALES 

The  characteristic  frequencies  in  the  general  dynamic  performance  of  ball 
bearings  normally  cover  a  wide  spectrum.  Frequencies  associated  with  the 
elastic  contact  phenomenon  are  on  the  high  end  while  the  shaft  speed  is  on 
the  low  end  of  the  spectrum.  Fortunately  for  most  operating  environments 
the  high  frequencies  are  several  orders  of  magnitude  greater  than  the  fre¬ 
quencies  of  the  order  of  shaft  revolution  speed  and  therefore  it  is  rela¬ 
tively  easy  to  constrain  the  motion  to  eliminate  all  the  high  frequencies 
in  order  to  look  at  the  low  frequency  pnenomenon  in  depth.  To  understand 
the  physical  mechanism  behind  each  frequency  the  basic  model  for  each  inter¬ 
action  and  the  fundamental  kinematics  of  ball  motion  must  be  reviewed. 

2.1  Ball/Race  Interaction 

The  normal  contact  loads  at  the  ball/race  interaction  are  computed  by  the 
Hertzian  elastic  contact  theory.  For  small  changes  in  load  the  non-linear 
Hertzian  load-deflection  relation  can  be  suitably  linearized. 

Now  if  the  ball  is  free  to  translate  arbitrarily  a  vibration  frequency 
corresponding  to  this  spring  will  be  observed.  This  has  been  discussed  in 
depth  by  Gupta  et  al  [6],  To  review  the  phenomenon  briefly  Figure  2-1  shows 
typical  axial  and  radial  components  of  ball  mass  center  as  it  travels  in 
its  orbit.  The  high  frequency  component  corresponds  to  the  elastic  contact 
phenomenon.  For  most  ball  bearings  this  frequency  can  be  in  the  range  of 
1  to  50  kHz.  A  large  bearing  will  normally  have  a  lower  elastic  contact 
frequency  due  to  heavier  balls. 

Figure  2-1  also  shows  a  distinct  low  frequency  whose  axial  and  radial 
components  are  180°  out  of  phase.  This  frequency  has  been  termed  as  a  kine¬ 
matic  frequency  [6]  since  it  is  a  strong  function  of  race  curvatures  and 
it  is  proportional  to  square  root  of  ball /race  contact  load.  Again  the 
actual  frequency  for  a  given  bearing  depends  on  the  geometry  of  bearing 
and  the  applied  conditions  but  the  general  range  is  from  about  500  Hz  to 
10  kHz. 

It  is  interesting  to  review  the  load  dependence  of  both  the  above  frequencies. 
The  kinematic  frequency,  by  definition,  has  a  square  root  variation  as  dis- 
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Fig.  2-1  Characteristic  Ball  Mass  Center  Vibration  Pattern 
In  an  Angular  Contact  Ball  Bearing 
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cussed  above  and  the  load  dependence  of  the  elastic  contact  frequency  can 
be  easily  shown  as  1/6  power  os  load: 


Hertz  contact  load 

Stiffness  k  ■  ■ 

Natural  frequency 


Q 


K63^2 


2/6 


Figure  2-2  shows  such  a  dependence  for  a  ball  bearing  with  12.7  mm  diameter 
balls  and  70  mm  pitch  diameter.  It  is  seen  that  as  the  contact  load  increases 
the  two  characteristic  frequencies  come  close  together.  Thus  either  at  high 
applied  loads  or  at  high  speeds,  the  distinction  between  the  two  frequencies 
may  not  be  as  clearly  seen  in  the  ball  acceleration  plots  as  it  is  seen  in 
Figure  2-1. 

2.2  Ball/ Cage  and  Race/ Cage  Interactions 

The  dynamics  of  ball/cage  interaction  can  be  quite  complicated  depending  on 
the  operating  conditions  and  the  ball  pocket  clearances.  For  a  simple  thrust 
loaded  condition  with  gravity  acting  along  the  bearing  axis,  the  cage  weight 
will  be  supported  uniformly  by  all  the  balls  and  therefore  ball/cage  contact 
will  be  established  in  each  pocket.  For  such  a  condition  the  elastic  contact 
spring  modeling  the  normal  contact  will  produce  a  vibratory  motion  if  the 
cage  is  allowed  to  translate  axially  relative  to  the  balls.  Such  a  phenomenon 
is  demonstrated  in  Figure  2-3.  Again  the  frequencies  will  be  generally 
high  compared  to  the  shaft  rotational  speed. 

In  most  instances  the  ball/cage  contacts  may  be  of  impulsive  nature  and  the 
total  time  of  contact  may  be  very  small.  A  typical  example  is  presented  in 
Figure  2-4.  Such  impulsive  contacts  will  produce  definite  discontinuities 
in  the  motions  of  both  the  balls  and  the  cage  and  considerable  care  will 
be  required  in  integrating  the  equations  of  motion. 

2.3  Low  Frequency  Components 

This  class  of  motion  basically  contains  frequencies  of  the  order  of  race 
angular  velocities.  The  components  of  interest  generally  include  ball 
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Fig.  2-2  Variations  of  computed  elastic  contact  frequency, 
oj  and  bearing  kinetic  frequency,  as  a  function 
o?  ball-race  contact  load 
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Fig.  2-3  Typical  Ball/Cage  Contact  Vibration 
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Fig.  2-4  Typical  Ball/Cage  Collisions  in  an  Angular  Contact  Ball  Bearing  with 
Combined  Axial  and  Radial  Load 


angular  velocity,  ball  orbital  velocity,  ball  pass  frequency  (rate  at  which 
a  given  point  on  the  race  is  passed  by  the  balls)  etc.  All  of  these  frequencies 
may  not  be  constant  depending  on  Che  operating  conditions.  In  fact,  the 
variation  in  ball  angular  velocity  may  sometimes  be  quite  large.  This  can 
be  easily  illustrated  if  the  equations  of  ball  angular  motion  are  written 
in  a  coordinate  frame  rotating  with  a  velocity  0  which  is  also  the  ball 
orbital  velocity: 


Iu),  -  lu)^0  m 
Io>2  +  Iu>2®  *  G^ 


Where  I  is  the  ball  moment  of  inertia,  u  is  the  angular  velocity  and  G 
is  the  applied  moment  vector. 


If,  for  instance,  the  applied  moment  is  zero,  then  will  be  constant  and 
0)2  and  will  be  governed  by 


o) 2  -  0)^9  "  0 
o>3  +  o>20  *  0 


or  u)„ 


u>3  9 


0  assuming  0  to  be  constant  and  also 


o)20 


0. 


Thus  both  o>2  and  will  have  a  cyclic  variation  with  a  frequency  of  0. 

The  nature  of  the  applied  moment  G  and  its  dependeance  on  ball  angular 
velocity  will  of  course  greatly  influence  the  ultimate  motion  of  the  balls. 


As  will  be  discussed  in  the  next  section,  a  knowledge  of  all  of  the  above 
characteristic  motions  is  necessary  if  any  constraints  on  the  ball  motion 
are  considered.  In  particular  the  knowledge  of  differences  between  the 
various  components  will  be  very  useful  and  two  classes  of  problems  may  be 
considered : 

Class  I:  All  frequencies  are  of  the  same  general  order 

Class  II:  There  is  a  distinct  difference  between  the  "very  low" 

and  "very  high"  frequency  components. 
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Class  I  will  be  relevant  for  a  relatively  large  bearing  operating  at  very 
high  speeds  and  Class  II  will  apply  to  most  small  to  moderate  size  bearings 
in  a  wide  operating  environment  and  large  bearings  operating  at  very  low 
speeds,  such  as  the  DMA  bearings.  For  the  Class  I  system  a  generalized 
motion  of  each  element  must  be  determined  while  certain  constraints  can  be 
imposed  on  Class  II  systems  depending  on  the  frequency  range  of  interest. 
Also,  some  constraints  may  apply  to  class  I  system  if  the  amplitudes  of 
motion  of  certain  components  are  very  small  compared  to  other  components. 

A  more  extensive  discussion  of  this  subject  is  presented  in  the  next  section 
in  terms  of  the  modifications  to  the  computer  program  DREB. 
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3.0  DREB  MODIFICATIONS  AND  RAPIDREB  DEVELOPMENTS 

The  available  Dynamics  of  Rolling  Element  Bearings  (DREB)  computer  program 
has  been  enhanced  for  several  added  capabilities  for  ball  bearing  performance 
simulation.  The  roller  bearing  capacilities  of  DREB  have  been  suppressed 
in  this  enhanced  version  named  RAPIDREB.  All  program  modifications  are 
centered  around  speeding  up  the  computation  and  examining  the  low  fre¬ 
quency  behavior  in  the  generalized  performance  simulation.  Hence  the  modi¬ 
fications  fall  in  two  catagories: 

1.  Suppression  of  High  Frequency  Content 

2.  Improved  integrating  algorithms. 

Aside  from  the  above,  a  minor  modification  to  extend  the  program  capabilities 
to  treat  ball  guided  cages  for  DMA  bearings  has  also  been  made.  This  section 
of  the  report  briefly  discusses  all  these  program  enhancements. 

3.1  Suppression  of  High  Frequency  Content 

As  discussed  in  the  preceding  section  there  are  certain  high  frequency 
characteristics  present  in  the  general  performance  simulation  of  ball  bear¬ 
ings.  A  careful  suppression  of  these  frequencies  is  necessary  to  efficiently 
simulate  the  low  frequency  phenomenon  specially  because  the  time  step  size 
in  the  numerical  integration  algorithm  is  determined  by  the  highest  fre¬ 
quency  present  in  the  system.  DREB  is  modified  in  two  ways  to  suppress  the 
high  frequency  motion,  e.  g.,  ball  motion  constraints  and  suitable  damping. 

3.1.1  Ball  Motion  Constraints 

The  very  high  frequency  motion  results  from  ball  vibration  between  the  two 
races.  Both  the  elastic  contact  frequency  and  the  kinematic  frequency 
represent  the  motion  of  ball  mass  center  relative  to  the  supporting  races, 
as  discussed  in  the  preceding  section.  If  the  applied  loads  on  the  bearing 
are  free  of  any  high  frequency  vibrational  loads,  the  amplitudes  of  the 
ball  vibration  are  generally  small  and  the  resulting  changes  in  contact 
loads  are  fairly  insignificant  when  compared  with  the  nominal  loads.  Thus 
it  may  be  reasonable  to  constrain  the  ball  such  that  the  ball  center  follows 
a  definite  path  which  satisfies  equilibrium  of  all  normal  contact  forces 
acting  on  the  ball  Including  the  centrifugal  forces.  The  traction  forces. 


being  quite  small  compared  to  the  normal  forces,  may  be  neglected  in  the 
equilibrium  equations.  Also  an  equilibrium  of  forces  is  carried  out  on  the 
inner  race  and  its  motion  is  constrained  to  satisfy  such  an  equilibrium. 

Such  constraints  are  reasonable  for  all  static  applied  loads  and  constant 
races  speeds.  They  are  also  applicable  to  time  varying  loads  and  speeds 
provided  the  frequency  of  such  variations  is  very  small  compared  to  the 
characteristic  ball/race  vibration  frequencies.  Thus  moderate  race  accelera¬ 
tions  and  synchronous  loading  due  to  unbalance  can  be  satisfactorily  treated 
within  the  realm  of  such  equilibrium  constraints. 

It  is  true  that  the  process  of  imposing  the  above  constraints  is  similar 
to  obtaining  the  conventional  quasi-static  solutions  at  each  time  step  and 
hence  a  set  of  non-linear  algebraic  equations  has  to  be  solved.  This  may 
appear  to  be  quite  time  consuming  but  even  for  relatively  large  steps, 
allowable  by  ball/race  traction  phenomenon,  the  changes  in  the  equilibrium 
conditions  are  rather  small  and  therefore  the  convergence  of  the  equilibrium 
equations  is  very  rapid.  Thus  in  spite  of  the  additional  effort  of  solving 
the  equilibrium  equations  at  each  step,  the  equilibrium  constraints  result 
in  a  substantial  reduction  in  the  overall  computational  effort,  primarily 
due  to  the  relatively  large  time  step  size  permissible  by  suppressing  the 
very  high  frequency  motion. 

3.1.2  Damping  Considerations 

There  may  be  some  damping  considerations  which  can  assist  in  eliminating 
certain  high  frequencies.  It  is  true  that  actual  damping  at  any  contact 
interface  in  the  bearing  is  quite  small  and  therefore  the  introduction  of 
any  damping  will  be  "fictitious"  with  the  prime  objective  of  eliminating 
the  high  frequency  response  without  significantly  altering  the  very  low 
frequency  behavior. 

Any  contact  interface  can  be  modeled  in  terms  of  a  spring  with  stiffness 
k  (N/M)  and  a  dashpot  with  damping  C  (NS/M)  as  shown  in  Figure  3-1.  The 
contact  stiffness  will  be  known  from  the  Hertzian  contact  parameters  at 
the  particular  contacts  and  the  damping  coefficient  will  be  the  required 
input.  Perhaps  specification  of  damping  will  be  best  in  terms  of  a  damping 
ratio  C/Cc>  where  Cc  is  called  critical  damping  [7]  corresponding  to  the 
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undamped  natural  frequency  of  the  system. 


Damping  ratio  ;  «  C/C 

c 

Where  C  *  2mw 
c  n 


Thus  c  can  be  specified  for  the  ball/race,  ball/cage  and  cage/race  inter¬ 
faces.  For  the  purpose  of  "real"  damping  5  will  be  very  small  for  most 
steels  but  for  certain  nonmetals  used  for  fabricating  the  cages  it  may  be 
appreciable.  In  any  event  any  non-2ero  value  of  c  will  only  effect  the 
general  response  in  the  neighborhood  of  the  undamped  natural  frequency 
[7]  and  the  response  at  very  low  frequencies  will  be  practically  unchanged. 

Each  contact  interface  in  the  bearing  is  modeled  in  terms  of  the  above  simpli¬ 
fied  vibration  model  and  the  computer  program  DREB  is  modified  accordingly. 

3.2  Integrating  Algorithms 

Algorithms  used  to  integrate  a  system  of  first  order  differential  equations 
may  be  either  of  implicit  or  explicit  type.  The  implicit  algorithms  are 
better  known  as  the  predictor  corrector  formulae  and  such  schemes  can  be 
very  efficient  if  no  discontinuities  and  high  frequency  oscillatory  motions 
are  present.  These  methods  make  use  of  a  set  of  existing  solutions  to 
extrapolate  to  an  expected  solution  at  the  current  time  step  and  this  is 
called  the  prediction  step.  With  the  predicted  values  the  derivatives  are 
computed  and  corrections  to  the  predicted  values  are  applied  in  an  iterative 
manner.  Normally  two  to  three  correction  iterations  are  enough  for  most 
non-stiff  systems.  If  no  high  frequency  motion  exists  in  the  bearing 
and  if  there  are  no  discontinuities  in  the  solutions  the  predictor-corrector 
process  can  provide  a  relatively  large  time  step  and  therefore  speed  up  the 
computation  in  the  desired  time  domain.  However,  if  discontinuities  are 
present,  such  as  ball/cage  or  race/cage  collisions,  the  extrapolation  process 
may  not  be  realistic  and  the  integration  process  as  a  whole  will  exhibit 
stability  problems.  For  such  cases  an  explicit  formula,  better  known  as  the 
Runga-Kutta  method,  which  makes  use  of  only  one  solution  to  march  to  the 


next  step  will  probably  be  the  best  choice. 

DREB  is  modified  to  incorporate  both  implicit  and  explicit  methods,  with 
appropriate  step  changing  criterion.  The  step  changing  process  generally 
depends  on  the  local  truncation  error  at  each  step.  An  estimate  of  this 
error  is  quite  straight  forward  for  an  implicit  method  but  it  may  be  quite 
difficult  for  an  explicit  method.  On  the  other  hand  treatment  of  variable 
size  steps  is  trivial  for  an  explicit  scheme  while  considerable  complications 
arise  in  the  extrapolation  process  of  the  predictor.  A  special  predictor- 
corrector  formula  in  terms  of  a  divided  difference  notation  is  therefore 
developed  to  efficiently  carry  out  the  extrapolation.  Some  suitable 
explicit  Runga-Kutta  methods  are  then  reviewed  with  particular  emphasis  on 
error  estimates  and  step  changing  rationale. 

3.2.1  Adams  Type  Predictor-Corrector  in  Divided  Differences 

A  predictor-corrector  algorithm  for  a  differential  equation  with  one  dependent 
variable  is  developed  in  this  section.  The  extension  of  the  algorithm  to  a 
system  is  quite  obvious.  Let  the  governing  differential  equation  be  written  as 

U  -  f(x,y),  (1) 

with  initial  condition  y0  ■  y(xQ). 

It  will  be  further  assumed  that  solutions  in  the  domain  0  <  x  <  x  have  been 

—  —  n 

obtained  by  some  explicit  algorithm  and  the  predictor-corrector  formula,  to  be 
developed  here,  will  advance  the  solution  from  xn  to  xn+^.  The  process  can  be 
written  as:  xn+l 

yn+l  "  yn  +  J  f(x*y)  dx  (2) 

y? 
n 

With  certain  polynomial  approximations  to  f(x,y),  the  required  algorithm  can 
be  developed.  The  formulation  presented  below  is  very  similar  to  the 
conventional  Adams-Bashford-Moulton  formulae. 

3. 2. 1.1  Predictor  Formula 

Let  the  solutions  (y^,  yn_^»  . yn-N^  and  ttie  derivatives  (^n*  fn-l ’  •  -fn_N) 

be  available  at  the  (N+l)  points  (x  ,  x  . . .  x  )  and  for  any  x,  let  the 

n  n-i  n-N 
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derivative  f(x)  be  a  Nth  order  polynomial  passing  through  the  points 


(V  xn-l . Vn5,  i,e': 


(x“xn-i+1>  f  [V  Xn-1' ' 


f(x)  -  f(xn)  +  £  (x-xn)  (x-xn _x)  . (x-xn_1+1)  f 

.4 

+  <x-xn)  (x-xn-1)  ....(x-xn_N)  f  [xn>xn_r...Xn^,C] 


"xn-i] 


Where  x  <  £  <  x  -  and  the  conventional  divided  difference  notion 
n-N  —  —  n+1 

[8]  is  used.  The  divided  difference  notation  will  be  convenient  for 

step  size  change  and  will  eliminate  the  need  for  interpolation  to  points  spaced 

uniformly  with  a  constant  grid  spacing. 

The  predictor  formula  is  obtained  by  substituting  the  polynomial  approximation, 
equation  (3),  in  equation  (2).  Thus 

N  i  j+1  j+1 

W  ■  +  A  £tV  Vl-'Vi1  \  *n  -  AJ 

i-0  J-0  (j+1) 

N+1  j+1  j+1 

+  ftxn’xn-l - xn-N,Cj  xn+l  ~  Xn  .  A.  (4) 

3  (j+1)  J 


Where  A^  are  the  coefficients  of  the  polynomial 

(x-xn)  (x-x^)  ••••(x-xn_i)  =  Ao  +  A1x  +  A^2  +  ....  Ai+1  xd 


Note  that  estimate  of  £  will  be  required  in  order  to  estimate  the  trunction 
error  term  of  equation  (4).  It  will  be  reasonable  to  assume  £  =  and  the 

corresponding  solution  as  the  predicted  value. 

From  computational  standpoint,  it  may  be  noted  that  as  the  order  N  increases, 
the  term  multiplying  the  polynomial  coefficient  A.,  may  experience  a  round 
off  problem,  which  can  be  eliminated  by  scaling  the  independent  variable  as 

*  Note  the  change  in  notation:  the  derivative  is  actually  a  function  of  y 
also,  notation  used  here  is  f(x)  ■  f(x,y). 
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X  .  ,  -X 
P+1  p 


so  that  the  polynomial 

(x-xn)  (x-xn  l) - (x-x^)^  Ap  (s"x*)  (S-Xn-1) - (s_xn-j) 

A 

where  x.  *  (x.  -  x  )/A 
J  j  P  P 

Also  if  p«*n,  the  integral  of  equation  (2)  takes  the  convenient  form: 


/* n+i  i+2  r1  *  * 

y{(x-xri)....(x-xn_j)}  dx  «  An  J  {(S-xn)  (S-xn-1) 


.  (s-xn-j)}  dS 


j+2  3+1  * 

-  A  I  ^k_ 

n  k-0  k+1 


fg 

Where  are  coefficients  of  the  polynomial 


*  cn“j+1 


(S_xn^ . (S-Xn-j^  =  A0  +  ^  S  +  - "An-j+l  S 


Substituting  (5)  and  (6)  in  equation  (4)  gives  the  final  predictor  formula 

N  £ 4.1  a  * 

y_ . ,  ■  y_  +  Z  f  [ x  ,  x  . , . . . .  x  ,  ]  A  Z  j 

n+1  n  i=»0  n  n_1  n_i  n  j-0  j+1 

N+2  N+1  * 

+  f[x  ,x  , ,  x  .,x  :  .]  A  Z  j  Cl 

n  n-1  n-N  n+1  n  j£j- 

* 

It  can  easily  be  seen  from  equation  (5)  that  the  coefficients  A^  will  only 
have  to  be  computed  if  the  grid  in  x  is  non-uniform. 

3. 2. 1.2  Corector  Formula 

For  the  general  corrector  formula,  let  f(x)  be  a  polynomial  of  degree  M, 

passing  through  the  (Mfl)  points  (xn+1>x  .x^. . .  •xn_M+1) .  Again  using  the 

divided  difference  notations 

M-l 

f(x)  =  f(xn+1)  +  Z  <*-xn+1)  (x-xn)  (x-x^-.-.Cx-x^^)  f[xn+1,xn,xn 

- xn-i 

+  (x-xn+l)  (x-xn) . (x-x^^)  ftxn+1,xn....xn_Wfl,  C] 
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Substitution  of  (8)  in  equation  (2)  gives  the  required  corrector  formula. 


M-l  i+1  i+1  1+1 

y  *  y  +  l  f  [x  .  ,,x . .  ,  ]  Z  n+l-n  B. 

n+1  n  i.-i  n+1  n  n’1  j-0  13+1 ) —  i 

+  f[xn+l’V  *  *  _ Ull  »  5  ]  1  4S,  -  *L+l 


n-M41 ’ 


j-o 


n+1  -n  B 


(j+1) 


j 


(9) 


Where  are  coefficients  of  the  polynomial 

(x-Xn+i)  (x-xn)  (xn_x)  ....(*-Vl)  =  B0  +  B1  x+  Bi+2  * 


i+2 


Again  an  estimate  of  the  truncation  error  will  require  the  value  of  £  which 
may  be  assumed  as  x^  ,  which  demands  that  one  additional  data  point  is 
available.  This  demand  is  easily  met  if  M  £  N. 


Scaling  the  independent  variable  again,  as  done  for  the  predictor  formula 
will  give  the  final  corrector  algorithm : 


yn+l 


M-l 

y  +  1  fix 

i*-l 


,x  , . . . . x  , ]  A 
n+1’  n’  n-i  n 


i+2 


+  f  [x 


n+1  *  n 


,x- - Xn-M+l,Xn-M^  A 


i+1  * 

^  !i 

1*0  j+l 
M+2  m+1 

I  B  * 

n  _i 


j-0 


(10) 


j+l 


As  for  the  predictor  formula,  the  coefficients  B^  will  only  have  to  be 
computed  for  non-uniform  grid.  Also  it  can  be  easily  seen  that  the  truncation 

error  will  employ  the  last  term  in  the  divided  difference  table  and  it  will  be 

.  ,  „  Mt-2 

proportional  to 


3. 2. 1.3  Computation  of  Coefficients 

In  order  to  implement  equations  (7)  and  (10)  a  formula  for  computing  the 

coefficients  of  the  polynomial 

P  =  (x— x  )  (x-x  .).... (x— x  ) 

n  n-i  n-p 

will  be  required. 


By  expanding  the  polynomial  for  p  ■  1,2,3, U  and  by  noting  the  algebraic 
similarities  in  the  expansions,  the  polynomial  can  be  written  in  the  following 
general  form: 


n+1  11+1  rs  M  V  P~k+2  p-k+3 

p  5  X^1  +  r  (-i)k  I  ,n  +1  J  X  (  X1 

k-l  i1-i  "_11  1  ij-ij+l  n'i2'fl 

P+1 

,Z  ..  xn-i, +1 
±k“1k-l+1  k 


p-k+j+1 

*  *  •  £  x  .  , ,  . 
i  4.1  n  — i  4*1 

rVi1  j 

(.id 


The  above  formula  can  be  coded  once  Che  maximum  value  of  p  Is  assigned. 


3. 2. 1.4  Step  Changing  Criterion 

The  step  size  can  be  changed  at  any  time  in  the  above  formulation  by  knowing 
the  truncation  error  at  the  preceding  step  size  and  the  allowable  limit.  Let 
us  assume  that  the  order  of  corrector  is  equal  to  that  of  the  predictor  (f*»N) 
and  hence  the  above  algorithm  is  a  (N+1)  step  process.  The  truncation  error 
of  the  predictor -corrector  step  will  now  be  given  by  the  truncation  error  of 
the  corrector  step  only  [9].  Thus  error  committed  in  arriving  at  xn+1>  frou 
xn  is  given  by: 


n+1 


.N+2 

4n 


fix. 


n+1’ 


•Xn-N1 


N+1 

I 


j-0 


(12) 


from  x  to  x  can  be  obtained  from  the 
n+1  n+2 

(13) 

Where  c  is  the  allowable  truncation  limit  and  C  (>  1)  is  a  constant,  which 

arbitrarily  reduces  the  computed  step  size  with  the  expectation  that  the  tolerance 

criterion  will  be  met  at  x  .  The  value  of  C  may  be  selected  primarily  by  ex- 

n+z 

perience,  it  may  vary  from  1.2  to  2.0.  Also  it  should  be  noted  that  any  time 
the  step  size  is  changed  all  coefficients  will  have  to  be  computed  for  N+1  steps 
(assuming  that  no  step  change  occurs  within  these  steps)  and  hence  this  added 
work  may  outweigh  the  savings  in  increased  step  size.  It  is  therefore  recommended 
that  the  step  size  be  changed  only  if  the  following  conditions  is  met: 


An  estimate  of  step  In  going 

relat ion : 


n+i  -rr 


n 

C 


1 

|N+2 
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>  a 


(14) 


Where  a  Is  a  constant,  again  selected  by  experience;  it  may  be  assumed  as  1.25 
as  a  first  approximation. 

Equation  (13)  can  also  be  used  to  reduce  the  step  size  at  any  time,  if  the 
truncation  limit  is  not  met. 

From  a  practical  standpoint,  upper  and  lower  limits  to  the  step  size  are  assigned 
and  the  objective  of  the  step  changing  rationale  is  to  keep  the  step  size  to  a 
maximum  within  this  range.  Also  if  the  truncation  limit  is  not  met  for  a  speci¬ 
fied  number  of  steps  then  further  integration  should  be  terminated. 

3. 1.2. 5  Change  of  Order 

A  criterion  for  change  of  order  is  not  very  clear  and  generally  a  trial  and 
error  approach  is  used.  At  a  given  time,  solutions  with  the  order  reduced  and 
increased  by  one  are  computed  and  the  order  giving  minimum  truncation  is  selected. 
Similar  to  the  step  size  the  maximum  and  minimum  bounds  on  the  order  should  also 
be  assigned.  The  maximum  order  can  generally  be  determined  by  examining  the 
divided  differences  which  will  generally  decrease  with  increase  of  order. 

However,  the  computed  values  may  show  an  increasing  trend  at  very  high  order, 
which  will  commonly  be  due  to  machine  round  off  errors.  Thus  the  maximum  order 
should  be  limited  by  this  phenomenon.  Selection  of  minimum  order  will  pri¬ 
marily  be  by  experience  with  the  problem  being  investigated. 

3. 1.2. 6  Computational  Considerations 

In  developing  a  computer  code  for  implementing  the  above  formulation,  certain 
considerations  will  be  helpful: 

1.  Set  a  maximum  limit  on  the  order  and  hence  the  number  of  solutions  to 
be  stored. 

2.  Compute  divided  differences  of  derivatives  at  the  first  call  and  upgrade 
table  at  subsequent  calls  by  adding  the  latest  solution  and  deleting  the 
earliest  solution. 
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3.  Check  grid  uniformity  and  compute  coefficients  (in  subsequent  calls) 
only  in  case  of  non-uniformity  resulting  from  step  size  change. 

U.  Check  maximum  allowable  order  by  examining  the  divided  difference  table 
and  make  sure  that  the  current  order  is  less  than  this  maximum. 

5.  Store  the  entire  divided  difference  table  in  disc  file  at  the  end  of  the 
run  for  subsequent  restarts  of  the  program. 

3.2.2  Explicit  Runga-Kutta  Tvoe  Algorithms 


As  discussed  earlier  an  explicit  formula  is  required  for  two  purposes. 

First  for  integrating  the  equations  in  the  domain  of  time  where  there  are 
either  discontinuities  or  large  cyclic  variations  at  relatively  high  fre¬ 
quency.  Secondly  an  explicit  formula  is  required  to  start  the  implicit 
predictor-corrector  process.  Also  it  is  required  that  some  estimate  of  the 
local  truncation  error  be  available  in  order  to  establish  some  criterion 
for  changing  the  step  size. 

The  existing  DREB  program  employs  a  Runga-Kutta-Merson  method  [9]  which 
does  estimate  the  truncation  error.  It  is  a  fourth  order  process  and 
requires  an  additional  function  evaluation  for  estimating  the  error.  Hence 
a  total  of  five  function  evaluations  per  step  are  used.  There  have  been 
some  questions  about  Che  accuracy  of  the  error  estimate  provided  by  this 
method  particularly  when  the  equations  are  non-linear.  Hence  methods  due 
to  Scranton  and  England  have  appeared  in  the  literature  [10] .  The  Scranton 
method  also  uses  five  function  evaluation  but  the  error  estimate  is  only 
realistic  for  a  single  equation  and  its  validity  to  a  system  of  equations 
is  questionable.  The  England  method  attempts  to  resolve  this  problem  by  an 
additional  function  evaluation  per  step  and  it  is  therefore  a  six  stage 
process.  Mere  recnetly  Fehlberg  [11]  has  provided  another  six  stage  formula 
which  has  been  of  considerable  use  in  recent  years.  Since  the  selection  of 
a  method  really  depends  on  the  problem  under  consideration,  which  in  the 
case  of  a  ball  bearing  can  be  quite  complicated  by  the  operating  conditions, 
both  the  formulae  due  to  England  and  Fehlberg  are  incorporated  in  the 
RAPIDREB  version  of  DREB. 
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All  of  the  explicit  formulae  can  be  best  summarized  in  terms  of 
tabular  notation  [12].  Thus  with  given  solution  yn  at  the  nth  step,  the 
solution  y  .  at  (n+l)th  step  is  given  by 


yn+l 


Yiki 


(15) 


where 


i-1 

ki  '  h  t<yn  +Z  1  *  1 . " 

J-l 


(16) 


and  the  estimate  of  local  truncation  error  is  given  by 


n+1 


I  Vk, 


i-1 


(17) 


The  values  of  S ,  Y  and  y*  for  the  three  methods  discussed  above  are  given  in 
Tables  3.1  to  3.3. 


All  of  the  above  methods  are  fourth  order.  For  improved  accuracy  a  higher 
order  method  may  sometimes  be  desirable,  e.g.,  in  a  detailed  simulation  of 
one  ball/cage  collision.  With  this  objective  a  fifth  order  Kutta-Nystrom 
method  [10]  and  a  sixth  order  (eight  stage)  Huta  method  [10]  are  also  incor¬ 
porated  in  the  computer  program.  The  coefficients  of  these  methods  are 
given  respectively  in  Tables  3.4  and  3.5.  For  these  methods  since  y*  is 
not  available  a  step  doubling  technique  is  used  to  compute  the  truncation 
error.  In  other  words  the  computation  at  each  step  is  done  twice,  first 
with  the  prescribed  step  size  and  then  with  half  the  prescribed  step  size. 
Thus  the  computation  with  these  methods (can  be  very  expensive  but  the  solu¬ 
tions  obtained  will  be  extremely  accurate.  Just  for  reference  purposes 
the  classical  fourth  order  Runga-Kutta  method,  the  fourth  order  Ralston- 
Runga-Kutta  method,  and  the  third  and  second  order  Runga-Kutta  methods  are 
also  incorporated  in  the  computer  program.  Again  a  step  doubling  technique 
is  used  in  all  these  methods  and  hence  considerable  care  must  be  exercised 
in  their  use.  The  coefficients  for  all  these  methods  are  fairly  well  known 
and  they  are  summarized  in  Tables  3.6  to  3.9  for  completeness. 
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TABLE  3.1 


COEFFICIENTS  FOR  THE  RUNGA-KUTTA-MERSON  METHOD 


TABLE  3.5 


COEFFICIENTS  FOR  THE SIXTH  ORDER  (EIGHT  STAGE)  HUTTA  METHOD 


TABLE  3.7 


COEFFICIENTS  OF  THE  FOURTH  ORDER  RALSTON-RUNGA-KUTTA  METHOD 

n  ■  4 


i 

6u 

Yi 

1 

0.17476028 

0 

0.40 

-  0.55148053 

3 

0.29697760 

0.15875966 

1.20553547 

4 

0.21810038 

-  3.05096470 

3.83286432 

0.17118478 

TABLE  3.8 


COEFFICIENTS  OF  THE  THIRD  ORDER  RUNGA-KUTTA  METHOD 

n  -  3 


«■!  pi  V  U.-. 


3.3  Modifications  for  Ball /Cage  Guidance  in  DMA  Bearings 
An  analytical  model  for  the  ball  guided  cages  has  been  presented  earlier 
[13].  The  model  assumes  a  conical  end  on  a  cylindrical  pocket  as  shown  in 
Figure  3-2.  It  is  really  a  purely  geometric  interaction  model  in  the  sense 
that  the  ball  center  is  located  with  respect  to  the  center  of  the  cone  base. 

-►  n 

This  is  shown  as  vector  r  in  the  pocket  frame  in  Figure  3-2. 

bp 


In  order  to  determine  the  magnitude  of 

it  will  be  convenient  to  transform  r£ 

bp 

(r,0,  z)  such  that 


interaction  with  the  conical  surface, 
into  cylindrical  coordinates 


'■V< 


>2  +  hi  >2 
bp2 


0  ■  arctan  (-r  P  /r  p  ) 

bp1  bp2 


(17) 


The  angle  0  measured  as  positive  rotation  about  the  z  axis  from  the  y  axis 
will  also  denote  the  contact  angle  for  the  cage  cone  contact.  It  should  be 
noted  that  if  the  cone  exists  on  the  inner  diameter  of  cage  the  positive  z 
axis  will  be  pointing  radially  inwards  and  the  y  axis  will  also  be  rotated 
by  180°. 


From  simple  geometry  it  may  be  shown  that  a  unit  vector  along  the  normal 
to  the  conical  surface  from  the  center  of  the  ball  is  given  by 


'bp 


cos<{i  sin  (9  +  0 ’ ) 
cos<t>  sin  (9  +  9’) 
sin$ 


(18) 


where  $  is  the  cone  angle  and  the  angle  0  is  measured  from  the  line  denoting  maxi¬ 
mum  interaction.  Thus  9’  -  0,  will  correspond  to  the  unit  vector  along  maximum 
interaction. 


The  interaction  6.  is  given  by 
bp 
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Fig.  3-2  Exaggerated  Pocket  Geometry  For  a  Ball 
Riding  Cage 


81 2825 
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(19) 


V0,)  “  (f  +  V  COS*  "  •  %l  -  f  (19) 

It  may  be  noted  that  since  e  ^  is  a  function  of  9'  the  interaction  will 

bp 

depend  on  9'and  it  will  be  necessary  to  determine  the  range  of  9'  for  which 

is  negative.  This  range  will  also  determine  the  length  of  interaction. 

Also  by  straightforward  geometry,  the  relative  position  of  the  ball  center 

and  the  length  of  the  conical  surface  will  determine  if  the  vector  e, 

bp 

falls  on  the  given  surface  or  not.  The  interaction  analysis  will  only  be 
performed  when  the  vector  falls  on  the  given  cage  surface.  Once  the 
point  of  contact  has  been  established  the  relative  slip  velocity  can  be 
readily  defined  in  terms  of  the  ball  and  cage  velocities. 


For  the  computed  geometric  interference,  the  normal  load  is  approximated  by 
either  a  Hertzian  point  or  line  contact,  depending  on  the  contact  arc 
length.  This  is  certainly  an  approximation  because  the  contact  between  a 
sphere  and  a  cone  is  not  really  Hertzian  and  the  true  solution  will  require 
a  generalized  three-dimensional  elastic  contact  analysis.  The  tractive 
forces  are  computed  by  substituting  the  calculated  slip  velocities  in  the 
prescribed  traction-slip  relation.  For  the  line  contact  the  contact  arc  is 
divided  into  several  incremental  areas  and  an  integration  is  performed  to 
compute  the  normal  and  traction  forces. 


3.4  Recommendations  for  Using  the  RAPIDREB  Version  of  DREB 


All  of  the  above  enhancements  are  incorporated  into  the  DREB  program  to 
create  the  new  version  RAPIDREB.  In  additio£,  with  the  objective  of  making 
the  program  efficient  for  ball  bearings,  the  roller  bearing  parts  and  some 
special  output  options  of  DREB  have  been  eliminated.  Depending  on  the  opera¬ 
tional  conditions,  the  computation  effort  in  RAPIDREB  may  be  significantly 
reduced  if  proper  options  are  exercised.  Some  useful  guidelines  to  this 
effect  are  therefore  presented  here: 


1.  If  elastohydrodynamic  lubrication  exists  it  is  recommended  that  a  few 
time  steps  be  run  on  the  original  DREB  program  and  typical  traction 
curves  be  printed  out.  Thus  a  ball  to  ball  variation  in  traction 
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may  be  studied  and  a  judgement  can  be  made  if  the  traction  behavior 
can  be  approximated  by  a  single  curve  and  thereby  eliminating  all  EHD 
computations  at  each  step.  If  so,  such  a  curve  must  be  approximated 
for  RAPIDREB. 

For  purely  thrust  loaded  bearings  the  ball  constraints  can  be  set 
such  that  the  load  computation  is  performed  only  once.  This  will 
significantly  effect  the  computational  effort. 

Bearings  with  a  combined  radial  and  thrust  load  may  also  employ 
ball/race  constraints  by  performing  an  equilibrium  analysis  at 
each  time  step.  This  may  result  in  an  increase  in  the  computational 
effort  per  step  but  the  total  effort  will  be  greatly  reduced  by 
permitting  the  maximum  step  size  to  be  much  larger  in  absence  of 
any  high  frequency  motions. 

Considerable  care  must  be  exercised  in  selecting  the  integration 
method.  By  extensive  running  of  RAPIDREB  it  is  found  tha*-  the 
predictor-corrector  scheme  performs  very  well  for  cageless  bear¬ 
ings  with  ball/race  constraints.  Hence  this  method  is  not 
recommended  for  most  bearings  with  cage,  where  large  discon¬ 
tinuities  in  the  solutions  may  be  present.  The  rehlberg 
and  England  formulae  do  sometimes  provide  an  increase  in 
permissible  step  size  but  an  increased  function  evaluation  at 
every  step  adds  to  the  effort  per  step.  The  benefit  of  these 
methods  may  only  be  investigated  for  prescribed  application  by 
a  few  trial  runs.  Any  of  the  higher  order  methods  provide  extreme 
accuracy  but  they  may  be  very  time  consuming  and  hence  these 
methods  should  only  be  used  in  very  small  time  domains. 

For  investigating  the  effects  of  small  perturbations  in  the  operat¬ 
ing  conditions,  a  steady  state  solution  under  nominal  conditions 
must  be  obtained  first  and  it  should  be  used  for  determining  the 
initial  conditions  for  subsequent  parametric  runs. 


6.  Proper  error  processing  control  must  be  used  to  process  time  limit 
errors  in  order  to  store  all  the  generated  data  files  in  the  event 
an  CP  time  limit  is  encountered.  This  is  normally  done  by  using 
the  EXIT  control  statement  on  the  CDC  systems. 

Appendix  A  contains  the  description  of  the  input  data  to  the  RAPIDREB  versions 
of  both  the  main  program  RDREB  and  the  plotting  program  RDREBP.  The  plot 
program  is  almost  identical  to  the  original  version  except  for  a  few  enhance¬ 
ments  for  handling  large  quantities  of  data  and  properly  truncating  the 
curves  if  they  exceed  the  plot  limits.  This  can  be  easily  understood  from 
the  input  description  in  Appendix  A. 


4.0  PERFORMANCE  SIMULATION  OF  DMA  BEARINGS 

Ball  bearings  used  in  DMA  (Despun  Mechanical  Assembly)  systems  are  generally 
very  large  bearings  operating  at  very  low  pure  thrust  loads  and  speeds. 
However,  they  are  designed  for  a  long  time  of  continuous  operation  (of  the 
order  of  ten  years).  The  lubricant  is  sealed  into  the  bearing  at  the  time 
of  assembly  and  hence  lubricant  degradation  or  starvation  is  not  an 
uncommon  problem.  Also  the  nature  of  friction  at  ball/cage  interface  has 
been  a  suspect  for  cage  instability  in  such  applications. 

Simulation  of  the  dynamic  performance  of  such  bearings  can  be  obtained  for 
extended  time  domains  by  using  the  equilibrium  constraints  on  ball  motion 
in  the  RAPIDREB  version  of  the  computer  program.  This  basically  eliminates 
the  ball/race  natural  vibration  and  thereby  allows  a  substantially  large 
time  step  to  integrate  the  differential  equations  of  motion  over  more  than 
a  shaft  revolution.  Such  simulations  are  the  subject  of  this  section. 

4.1  Bearing  Geometry 

Two  different  geometrical  configurations  of  the  DMA  bearings  are  investi¬ 
gated.  Figure  4-1  shows  the  geometrical  details  of  a  100  mm  bore  angular 
contact  bearing  with  the  cage  guided  on  the  inner  race  and  Figure  4-2  dis¬ 
plays  the  design  of  a  150  mm  bore  ball  bearing  with  the  cage  guided  on  the 
balls.  The  guidance  surface  is  of  a  conical  shape  and  it  is  attached  to 
the  cylindrical  pockets  on  the  inner  diameter  of  the  cage.  Details 
of  the  cage  geometry  are  shown  in  Figure  4-3. 


4.2  Operating  Conditions 

Both  bearings  operate  at  relatively  low  pure  thrust  loads  and  speeds, 
which  are  typical  of  actual  application.  The  operating  condition  assumed 
for  the  simulations  to  be  presented  here  are: 


100  mm  bearing : 


150  mm  bearing: 


Thrust  load  ** 

Inner  race  speed  * 

Thrust  load  ■ 

Inner  race  speed  * 


178  N 
60  RPM 

534  N 
62.5  RPM 
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Race  Guided  Cage 


The  outer  race  is  assumed  to  be  stationary  in  both  cases  and  the  operating 
temperature  of  both  bearings  is  assumed  to  be  15°C.  The  motion  of  the  cage 
is  considered  under  the  influence  of  gravity  with  the  bearing  axis  oriented 
both  along  and  normal  to  the  direction  of  gravity.  This  is  really  the 
only  factor  which  will  promote  cage  motion  and  its  interaction  with  the 
balls.  The  two  orientations  of  the  bearing  with  respect  to  the  direction 
of  gravity  do  represent  laboratory  conditions. 

4.3  Ball/Race  Traction  Models 

Traction  behavior  of  typical  lubricants  used  in  DMA  bearings  is  not  fully 
known  and  a  number  of  assumptions  in  modelling  the  behavior  are  generally 
made.  For  the  present  application  a  lubricant  called  Apiezon  C  +  1.5% 
antitoxidant  [14,  15]  is  considered  with  the  following  properties  under 
ambient  pressure  and  at  a  reference  temperature  of  15°C. 

Viscosity  =  0.38425  Pa.S 

-9 

Pressure-Viscosity  Coefficient  =  8.4234  x  10  1/Pa 

Temperature-Viscosity  Coefficient  *  6G00°K 

Thermal  Conductivity  *  0.10  N/S/°C 

Since  the  traction  parameters  are  not  known,  the  following  parameters,  which 
actually  correspond  to  a  MIL-L-7808  type  oil,  are  arbitrarily  assumed  for 
the  Type  I  traction  model  discussed  by  Gupta  et  al  [16]. 

Viscosity  Parameter,  u*  *  0.58889  Pa 

_9 

Pressure-Viscosity  Parameter,  a*  «  5.22136  x  10  1/Pa 

_2 

Temperature  Viscosity  Parameter,  8*  *  2.48550  x  10  1/°K 

Also  typical  of  the  DMA  bearing  environment  a  lubricant  starvation  factor 
(distance  of  lubricant  maniscus  from  the  edge  of  the  Hertzian  contact  zone 
divided  by  the  major  contact  half  width)  of  2  is  assumed. 

With  the  above  lubricant  parameters  and  the  operating  conditions  outlined 
earlier  the  computed  elastohydrodynamic  traction  coefficients  for  the  two 
bearings  are  fairly  linear  with  the  slide-to-roll  ratio,  in  the  range  of 
0  to  0.10,  as  shown  in  Figure  4-4.  The  lubricant  film  thicknesses  are 
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Fig.  4-2  Geometry  of  the  150  mm  Df-'A  Ball  Bearing  with  Ball  Guided  Cage 


also  listed  in  the  figure  and  from  these  numbers  it  may  be  concluded 
that  the  assumption  of  full  film  elastohydrodynamic  lubrication  even  under 
starved  conditions  may  be  fairly  reasonable.  If  the  data  of  Figure  4-4  is 
plotted  as  a  function  of  the  sliding  speed  or  slip  velocity,  then  all  of 
the  data  lies  on  one  line  with  a  slope  of  0.0410  S/M;  this  is  because  the 
roll  speed  is  different  for  each  line  in  Figure  4-4.  Thus  a  linear  traction 
curve  with  a  slope  of  0.0410  S/M  is  assumed  for  both  bearings. 

Friction  at  ball/cage  and  cage/race  contacts  is  another  unknown  and  for  the 
present  investigation  a  constant  friction  coefficient  of  0.010  is  arbitrarily 
assumed . 

4.4  Performance  Simulations  of  100  mm  DMA  Bearing 

The  first  simulation  with  the  RAPIDREB  program  is  obtained  for  the  100  mm, 
inner  race  guided  DMA  ball  bearing  operating  with  gravity  oriented  along 
the  bearing  axis.  In  such  a  orientation  the  cage  moves  axially  and  rests 
on  all  the  balls  to  balance  its  own  weight.  Typical  ball/cage  force  varia¬ 
tion  is  shown  in  Figure  4-5.  The  mean  force  is  approximately  equal  to  the 
cage  weight.  What  is  interesting  to  notice  is  the  settling  in  of  the  ball/cage 
contact  angle.  The  observed  distinct  frequency  in  this  pattern  is  about 
44.6  Hz  and  it  is  probably  related  to  some  kinematic  phenomenon  in  the  bearing; 
it  is  really  much  lower  than  the  ball/cage  contact  frequency  which  is  of 
the  order  of  1  kHz.  Such  a  cyclic  variation  also  gets  transmitted  to  the 
cage  angular  acceleration  about  the  bearing  axis  as  shown  in  Figure  4-6, 
and  hence  to  the  skid  parameter  as  shown  in  Figure  4-7.  The  parameter 
A 2  is  defined  as: 

X  ,  Cage  Angular  Velocity 

2  Race  Angular  Velocity 

From  Figure  4-7  it  is  also  seen  that  the  cage  does  not  experience  any  signifi¬ 
cant  whirl.  This  is  indicated  by  the  parameter  A^  in  Figure  4-7  which  is 
defined  as: 

^  m  Cage  Mass  Center  Whirl 

1  Cage  Angular  Velocity 
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Fig.  4-6  Dimensionaless  Cage  Angular  Accelerations  for  the 
100  mm  DMA  Bearing  with  Gravity  Acting  Along  the 
Bearing  Axis.  Scale  -  1.319  x  10'  RPM/S. 
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The  sudden  change  from  zero  whirl  to  some  vibratory  motion  indicates  cage 
radial  motion  beyond  a  certain  machine  round  off  threshold  set  in  RAPIDREB. 
In  reality  the  entire  whirl  pattern  will  look  like  what  is  seen  in  steady 
state  in  Figure  4-7,  but  the  magnitude  of  the  radial  motion  for  all  practi¬ 
cal  purposes  will  be  negligible  and  no  interaction  with  the  guiding  race 
will  be  seen.  The  bearing  torque  variations  are  seen  in  Figure  4-8.  Again 
the  cyclic  variation  corresponding  to  about  44.6  Hz  is  seen  in  both  the 
outer  and  inner  race  torque.  The  damping  of  the  cyclic  variation  is  due 
to  ball/race  and  ball/cage  friction.  Under  actual  laboratory  conditions 
such  cyclic  torque  variations  may  exist  in  steady  state  since  the  bearing 
is  constantly  subjected  to  subtle  disturbances  which  may  be  adequate  to 
continually  excite  the  cyclic  motion. 

The  ball  motion  indicates  that  the  conventional  race-control  hypotheses 
do  not  hold  when  the  ball/race  friction  is  quite  low  as  in  the  present  case. 
As  shown  in  Figure  4-9  the  angular  velocity  vector  tends  to  orient  itself 
along  the  bearing  axis  and  hence  the  x  component  continually  increases 
while  the  z  component  decreases.  Also  a  small  component  about  the  y  axis 
is  developed  indicating  gyroscopic  slip.  The  resulting  spin-to-roll  ratios 
are  shown  in  Figure  4-10.  Such  deviations  of  ball  velocity  from  the  race- 
control  hypothesis  have  also  been  observed  earlier  [17]. 

When  the  bearing  is  oriented  vertically  such  that  the  cage  weight  now  acts 
in  a  plane  perpendicular  to  the  bearing  axis,  the  cage  tends  to  move  radi¬ 
ally  and  this  causes  substantial  variations  in  ball/cage  interactions  as 
the  ball  travels  around  its  orbit.  Also  significant  interaction  with  the 
guiding  race  is  observed  and  some  coning  motion  (rotation  about  transverse 
axes  of  the  cage)  is  also  noted.  Figure  4-11  shows  the  ball/cage  force 
variation.  The  high  frequency  content  is  generally  the  ball/cage  elastic 
contact  frequency  and  the  low  frequency  variation  corresponds  to  the  ball 
orbital  velocity;  the  simulation  time  in  Figure  4-11  corresponds  to  about 
one  shaft  revolution  or  approximately  half  of  ball  orbit.  The  cage/race 
forces  are  shown  in  Figure  4-12.  The  steady  magnitude  of  the  force  is  of 
the  order  of  cage  weight  although  not  exactly  equal  to  it  since  the  ball/ 
cage  forces  also  enter  in  the  cage  equilibrium.  Although  not  very  clear 
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Fig.  4-8  Bearing  Torque  Variations  for  the  100  mm  DMA 
Bearing  with  Gravity  Acting  Along  the  Bearing 
Axis 
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Fig.  4-10  Ball/Race  Spin-to-Roll  Ratio  for  the  100  mm 
DMA  Bearing 
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Fig.  4-11  Ball/Cage  Force  Variation  for  the  100  mm 
Bearing  with  Gravity  Acting  Normal  to  the 
Bearing  Axis 
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from  Figure  4-12,  there  is  an  appreciable  difference  between  the  contact 
forces  at  the  two  land’,  which  is  due  to  the  coning  motion  of  the  cage. 

This  is  better  seen  in  Figure  4-13  where  the  cage  angular  accelerations 
are  plotted.  Most  of  this  motion  is  clearly  cyclic  and  it  has  both  the 
high  and  low  frequency  content  present  in  it.  The  high  frequency  corres¬ 
ponds  to  the  ball/cage  and  cage/race  elastic  contacts  while  the  low  fre¬ 
quency  is  similar  to  what  is  seen  earlier  in  Figure  4-5.  Both  the  y  and  z 
components  do  show  the  low  frequency  but  it  is  more  clear  in  Figure  4-14, 
where  the  angular  velocities  are  plotted.  The  magnitude  of  this  low  fre¬ 
quency  is  about  76.4  Hz  compared  to  the  frequency  of  44.6  Hz  shown  earlier 
with  the  bearing  more  symmetrically  oriented. 

In  spite  of  the  complicated  motion  the  cage  does  seem  to  be  fairly  stable 
with  practically  no  whirl  in  steady-state  as  indicated  by  A^  in  Figure  4-15. 

Also  seen  in  this  figure  is  the  skid  parameter  A^,  which  is  also  fairly 
stable.  The  bearing  torque  variations  are  shown  in  Figure  4-16.  These 
variations  do  contain  both  the  high  and  low  frequency  cyclic  variations , 
although  they  are  not  very  clear  in  the  figure.  The  ball  angular  velocities, 
as  seen  in  Figure  4-17,  are  similar  to  those  seen  earlier.  Also,  the  devia¬ 
tions  from  the  race-control  hypothesis  are  clear  in  terms  of  the  spin-to-roll 
ratios  plotted  in  Figure  4-18.  In  addition  to  this  it  may  be  interesting  to 
note  the  ball  orbital  acceleration  pattern  shown  in  Figure  4-19.  The  noise 
in  the  ball  orbital  acceleration  correlates  well  with  the  noise  observed  earlier 
in  the  skid  parameter  in  Figure  4-15.  Again  the  high  frequency  content 
comes  from  the  ball/cage  elastic  collisions,  which  are  now  continually  varying 
as  the  ball  travels  around  its  orbit. 

4.5  Performance  Simulations  of  the  150  mm  DMA  Bearing 

The  general  behavior  of  the  150  mm  DMA  bearing  with  ball  guided  cage  is 

relatively  more  noisy  than  the  100  mm  bearing.  With  the  gravity  oriented 

along  the  bearing  axis  the  ball/cage  pocket  force  variation  is  shown  in 

Figure  4-20.  Once  again  the  mean  force  balances  the  cage  weight  and  the 

ball/cage  contact  angle  demonstrate  some  cyclic  variations  as  is  seen  for 

the  100  mm  bearing  (Figure  4-5).  The  frequency  of  this  component 

is  Oow  37.4  Hz  which  is  not  very  different  from  44.6  Hz  observed  with  the  100  mm 

bearing.  The  resulting  cage  motion  is,  however,  substantially  noisier. 
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Fig.  4-13  Dimensionless  Angular  Accelerations  of  the 

100  mm  Bearing  Cage  When  Gravity  Acts  Normal 
to  the  Bearing  Axis.  Scale  *  1.319  x  10  RPM/S 
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Fig.  4-15  Cage  Whirl  and  Skid  Parameters  for  the  100  mm 
Bearing  with  Gravity  Acting  Normal  to  Bearing 
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Fig.  4-16  Variation  in  Bearing  Torque  in  Case  of  the 
100  mm  Bearing  with  Gravity  Acting  Normal 
to  the  Bearing  Axis 
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Fig,  4-20  Ball/Cage  Interaction  in  the  150  mm  Bearing  With 

Ball  Guided  Cage  Operating  Horizontal  with  Gravity 
Acting  Along  the  Bearing  Axis 


Figure  4-21  shows  cage  mass  center  accelerations.  It  is  clear  that  both 

radial  and  orbital  components  develop  beyond  the  threshold  of  machine  round 
off  errors.  Although  the  resulting  whirl  motion  is  substantial  in  terms  of 

angular  velocity,  the  whirl  radius  is  still  small  and  no  appreciable  inter¬ 
action  with  the  pocket  cone  takes  place.  The  cage  whirl  is  indicated  by  A^ 
variation  in  Figure  4-22.  It  is  easily  seen  that  once  the  cage  starts 
whirling  the  whirl  velocity  oscillates  by  large  magnitudes;  this  may  indi¬ 
cate  some  form  of  instability.  Cage  angular  velocity,  indicated  by  A^, 
is,  however,  fairly  steady  after  the  dampening  of  the  initial  cyclic  varia¬ 
tions.  The  bearing  torque  variations  are  very  similar  to  those  for  the  100 
bearing,  as  seen  in  Figure  4-23. 

The  general  motion  of  the  ball  is  almost  identical  to  that  observed  in  the 
100  mm  bearing.  The  ball  angular  velocity  vector,  see  Figure  4-24,  tends 
to  line  up  with  the  bearing  axis  and  the  race-control  hypotheses  do  not 
hold,  as  seen  by  the  spin-to-roll  ratios  in  Figure  4-25. 

When  the  150  mm  bearing  with  ball- guided  cage  is  run  with  gravity  oriented 
normal  to  the  bearing  axis,  significant  noise  in  the  performance  simulations 
is  produced  by  the  repeated  collisions  in  the  cage  pocket.  Figure  4-26 
shows  typical  collisions  with  the  cylindrical  part  of  the  pockets  and  Figure 
4-27  shows  interaction  at  the  conical  guidance  surface.  Such  repeated  colli¬ 
sions  produce  considerable  noise  in  both  the  ball  and  cage  motions.  The 
ball  angular  velocity  variations  are  shown  in  Figure  4-28  and  the  orbital 
accelerations  of  the  ball  from  the  repeated  impacts  is  shown  in  Figure  4-29. 
The  angular  velocity  of  the  cage  indicates  substantial  coning  motion  as 
seen  by  the  y  and  z  components  in  Figure  4-30.  The  noise  generated  from  the 
ball-cage  collisions  is  also  transmitted  to  the  cage  mass  center  accelera¬ 
tions  as  shown  in  Figure  4-31.  Excessive  cyclic  whirl  of  the  cage  is  pro¬ 
duced  although  no  significant  steady  whirl  is  observed  over  the  short  time 
of  this  performance  simulation.  This  is  seen  by  the  A^  variation  in 
Figure  4-32.  As  may  be  expected  the  skid  parameter,  A^,  does  pick  up  the 
influence  of  the  excessive  ball/cage  collisions. 
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Fig.  4-21  Dimensionless  Cage  Mass  Center  Accelerations 

Resulting  From  the  Noisy  Ball/Cage  Interaction 
for  the  150  mm  Bearing.  Scale  *  1.198  x  10  M/S 
1.03  x  107  RPM/S  cQ 
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Fig.  4-23  Bearing  Torque  Variation  for  the  150  mm  Bearing 
with  Gravity  Acting  Along  the  Bearing  Axis 
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Fig.  4-25  Ball/Race  Spin  for  the  150  am  DMA  Bearing 
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Fig.  4-29  Ball  Orbital  Accelerations  Resulting  From  the 
Excessive  Ball/Cage  Interaction  for  the  150  mm 
DMA  Bearing 
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Fig.  4-30  Dimensionless  Angular  Velocity  of  the  Ball  Guided 
Cage  in  the  150  mm  Bearing  with  Large  Ball/Cage 
Interaction.  Scale  ■  9.918  x  10J  RPM 
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Fig.  4-32  Cage  Whirl  and  Skid  Parameters  for  the  150  mm 
Bearing  with  Ball  Guided  Cage  and  with  Gravity 
Acting  Normal  to  the  Bearing  Axis 


T 


5.0  ENGINE  BEARING  PERFORMANCE  SIMULATIONS 

The  enhanced  capabilities  of  RAPIDREB  also  apply  to  the  high~load  and  high¬ 
speed  engine  bearings  by  providing  bearing  performance  simulation  over 
several  shaft  revolutions.  For  a  given  bearing  and  a  set  of  operating 
conditions  some  typical  performance  simulations  are  presented  in  this  section. 

5.1  Bearing  Geometry 

The  geometry  of  the  bearing  is  shown  in  Figure  5-1.  It  is  a  100  mm  bore 
ball  bearing  with  a  steel  cage  guided  on  the  inner  race.  The  design  corre¬ 
sponds  to  an  actual  test  bearing  at  the  Air  Force  Propulsion  Laboratories. 

5.2  Lubricant  Traction  Models 

Conventional  lubricant  with  the  MIL-L-7808  specification  is  assumed.  The 
traction  behavior  of  this  lubricant  has  been  studied  fairly  extensively 
and  the  model  in  the  original  DREB  program  has  been  demonstrated  to  show  a 
reasonable  fit  with  the  experimental  data.  The  model  is  also  built  into 
RAPIDREB  and  appropriate  input  option  is  exercised  to  select  the  MIL-L-7808 
model . 


For  the  ball/cage  and  cage/race  interfaces  a  hypothetical  traction  model 
(see  Appendix  C)  is  used  to  compute  the  friction  forces  when  a  metal  contact 
takes  place. 


5.3  Operating  Conditions 
To  simulate  typical  test  conditions 
assumed : 

Axial  load  - 
Radial  load  = 

Inner  race  speed  * 

Outer  race  speed  * 
Operating  temperature  = 


the  following  operating  conditions  are 

18,000  N 
4,500  N 
20,000  RPM 
0 

330°K 


Gravity  acts  normal  to  the  bearing  axis,  but  in  this  case  it  really  plays 
no  appreciable  role  because  the  applied  forces  are  quite  high.  Radial  and 
axial  equilibrium  contraints  on  ball  motion  are  assumed  to  eliminate  the 
high  frequency  ball/race  vibrations.  Also  the  quasi-static  solutions  are 
used  to  determine  the  initial  conditions. 
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5.4  Performance  Simulations 

Performance  simulations  for  the  engine  bearing  at  the  above  conditions  are 
obtained  over  seven  shaft  revolutions  and  hence  definite  steady  state 
behavior  can  be  easily  understood.  Typical  initial  parameters  are  illustrated 
in  the  computer  output  presented  in  Appendix  C. 

As  might  be  expected  in  case  of  combined  axial  and  radial  load,  the  ball 
load  and  contact  angles  go  through  a  cyclic  variation  with  a  frequency 
corresponding  to  the  ball  orbital  velocity.  This  is  seen  in  Figure  5-2, 
where  over  three  wavelengths  of  variations  corresponding  to  over  seven  shaft 
revolutions  are  shown.  It  is  also  seen  in  this  figure  that  the  race- 
control  hypotheses  do  not  hold  heie  also  and  definite  spin  velocities  on 
both  races  develop  in  steady-state.  Ball  angular  velocities  demonstrate 
an  interesting  effect,  see  Figure  5-3.  Since  the  initial  quasi-static 
solutions  do  not  allow  gyroscopic  slip  and  exact  equilibrium  of  moments  is  not 
satisfied,  the  ball  immediately  tends  to  slip  about  the  transverse  y 
axis  due  to  the  gyroscopic  moments.  This  alters  the  x  and  z  components 
also  and  it  takes  approximately  two  shaft  revolutions  for  the  angular 
velocities  to  develop  some  steady  state  pattern  and  satisfy  the  gyroscopic 
slip  as  permitted  by  the  lubricant  characteristics.  The  ball  cage  inter¬ 
actions  show  that  about  one  collision  in  each  pock  takes  place  per  revolu¬ 
tion  of  the  cage  and  the  ball  drives  the  cage  (ball/cage  contact  angle  of 
180°)  for  part  of  the  revolution  and  gets  driven  by  the  cage  (contact  angle 
of  zero)  for  the  remaining  part.  This  is  seen  in  Figure  5-4.  The  ball/ 
cage  forces  seen  in  this  figure  are  really  the  collision  forces  and  the 
hydrodynamic  forces,  being  very  small  compared  to  the  collision  forces, 
are  not  seen  in  the  figure.  However,  the  ball/cage  approach  relative  to 
the  radial  pocket  clearance  clearly  demonstrates  the  transition  from  hydro- 
dynamic  to  metal  contact. 

The  race/cage  interaction  is  very  dynamic  and  interesting  to  note.  Figure 
5-5  shows  the  race/cage  force  variations  at  the  two  lands  (labelled  1  and  2) 
on  either  side  of  the  balls.  Intially  both  lands  show  identical  collision 
forces  but  when  appreciable  coning  motion  of  the  cage  develops,  the  forces 
on  the  two  lands  begin  to  differ;  this  is  more  clear  from  the  traction 
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5-3  Dimensionless  Ball  Angular  Velocity  Variation  for  the 
100  mm  Engine  Bearing.  Scale  -  7.84  x  104  RPM 
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Repeated  Cage/Race  Collisions  in  the  100  mm  Engine  Bearing 


force  curves.  The  ultimate  behavior  to  note  is  that  the  race/cage 
collisions  are  really  very  close  together  in  time  in  steady  state.  Also 
the  magnitude  of  the  force  is  fairly  large  (*>  1800  N) .  This  means  that 
the  cage  will  steadily  be  in  contact  with  the  race  with  a  relatively  large 
force,  and  some  cage  wear  problems  may  be  expected.  This  is  seen  by  the 
steady  circular  orbit  of  the  cage  mass  center  in  Figure  5-6  with  the  orbit 
radius  equal  to  the  cage/race  radial  clearance.  This  effect  has  indeed 
been  confirmed  experimentally  [18]  and  the  experimental  results  will  soon 
be  documented  elsewhere  [19]. 

Corresponding  to  the  above  race/cage  collisions  considerable  radial  chatter 
of  the  cage  mass  center  is  observed,  as  shown  in  Figure  5-7  in  terms  of  the 
cage  mass  center  velocities.  Excessive  whirl  is  indicated  by  the  increased 
orbital  velocity.  The  coning  motion  is  demonstrated  by  the  angular  velocity 
variations  in  Figure  5-8.  The  deviation  of  cage  angular  velocity  from  the 
initial  epicyclic  angular  velocity  (see  component  x)  indicates  substantial 
skid  in  the  bearings  due  to  excessive  rubbing  of  the  cage  at  the  cage/race 
interrace.  Both  cage  whirl  and  skid  are  better  understood  in  terms  of  the 
parameters  X^  and  X2  respectively,  as  shown  in  Figure  5-9. 

Finally,  the  variations  of  bearing  torque,  power  loss,  and  the  cumulative 
load  slip  integral  are  shown  in  Figure  5-10.  The  initial  bump  in  the  torque 
curves  correspond  to  the  gyroscopic  slip  of  the  balls  as  discussed  above 
and  the  steady  noise  is  a  result  of  the  ball/cage  and  race/cage  collisions. 
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Fig.  5-9  Cage  Whirl  and  Skid  Parameters  for  the  100  mm  Engine 
‘  Bearing 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  existing  Dynamics  of  Rolling  Element  Bearings  (DREB)  computer  program 
has  been  enhanced  to  selectively  suppress  the  very  high  frequency  vibra¬ 
tory  motion  of  the  balls  in  an  angular  contact  ball  bearing  and  thereby 
provide  the  added  capability  of  investigating  the  low  frequency  phenomena 
in  some  depth.  With  such  an  enhancement  a  considerable  increase  in  the 
maximum  permissible  time  step  size  has  led  to  bearing  performance  simula¬ 
tion  over  several  shaft  revolutions  within  reasonable  computing  effort. 
The  updated  program  version  is  referred  to  as  RAPIDREB.  Some  of  the 
essential  RAPIDREB  capabilities  include: 

1.  Radial  and  axial  equilibrium  constraints  on  the  balls  to 
eliminate  the  high  frequency  ball/race  vibration. 

2.  Introduction  of  fictitious  damping  at  each  contact  interface 
to  selectively  damp  the  high  frequencies. 

3.  Incorporation  of  a  number  of  Runga-Kutta  type  explicit  formulae. 

4.  A  predictor-corrector  algorithm  for  conditions  under  which  no 
substantial  cyclic  components  and  discontinuities  in  the  solu¬ 
tions  exist. 

5.  Option  to  change  the  integrating  algorithm  at  any  point  in 
time. 

6.  Treatment  of  ball  guided  cages  for  very  low  speed  ball  bearings. 

7.  Built  in  options  for  synchronous  radial  loading  due  to  unbal¬ 
ance  or  any  rotating  loads. 

Any  of  the  above  options  can  be  selected  or  the  program  can  be  run  in  a 
completely  generalized  mode  identical  to  the  original  DREB  program.  It 
is  shown  that  depending  on  the  operating  conditions  substantial  computing 
effort  can  be  saved  by  properly  selecting  the  constraints  and  the  inte¬ 
grating  method.  By  exercising  RAPIDREB  over  a  number  of  trial  runs  it  is 
concluded  that  the  predictor-corrector  type  of  formula  is  only  suitable 
when  there  are  no  discontinuities  and  high  frequency  vibratory  motions  in 
the  system.  Thus  after  suppressing  the  very  high  frequencies,  the 


predictor-corrector  schemes  offer  the  benefit  of  very  large 
step  sizes,  only  for  the  cases  where  no  discontinuities  associated  with 
ball/cage  or  race/cage  collisions  are  present.  This  is  because  extrapola¬ 
tion  past  a  discontinuity  results  in  larger  errors  in  the  prediction  process. 
Thus  these  schemes  are  useful  only  for  bearings  either  without  cage  or 
with  basically  hydrodynamic  interaction  at  the  ball/cage  and  race/cage 
interfaces.  Considerable  care  in  selecting  the  integrating  method  is  therefore 
necessary. 

RAPIDREB  is  used  to  produce  performance  simulations  of  two  very  low  speed, 
lightly  loaded  DMA  (Despun  Mechanical  Assembly)  ball  bearings  and  a 
typical  high-load,  high-speed  ball  bearing  for  engine  applications.  The 
DMA  bearings  have  pure  thrust  load  and  therefore  gravity  is  the  only 
dominant  factor  as  far  as  the  cage  is  concerned.  Hence  both  bearings  are 
oriented  such  that  gravity  is  either  parallel  or  perpendicular  to  the 
bearing  axis.  The  first  of  the  DMA  bearings  is  a  100  mm  bore  ball  bearing 
with  inner  race  guided  cage.  The  steady  state  performance  simulation 
obtained  by  integrating  the  equations  of  motion  over  more  than  a  shaft 
revolution  indicate  that  when  gravity  acts  along  the  bearing  axis  the  cage 
contacts  all  the  balls  and  balances  its  own  weight  against  the  balls.  No 
substantial  cage  whirl  is  observed  in  the  simulation  and  the  bearing  opera¬ 
tion  is  generally  stable.  When  weight  of  the  cage  acts  normal  to  the  bear¬ 
ing  axis  variable  ball/cage  and  cage/race  interactions  are  observed  as  the 
balls  travel  in  their  orbit.  This  is  basically  due  to  the  eccentric  loca¬ 
tion  of  the  cage.  An  appreciable  coning  motion  of  the  cage  is  also  observed 
in  this  case  but  again  no  appreciable  cage  whirl  or  skid  is  observed  and  the 
operation  is  generally  stable.  The  second  of  the  DMA  bearings  is  a  150  mm 
bore  ball  bearing  with  a  ball  guided  cage  with  a  conical  guidance  surface 
attached  to  the  nominally  cylindrical  pocket.  Again  the  first  symmetric 
case  when  the  gravity  acts  along  the  bearing  axis  is  similar  to  the  100  mm 
bearing  and  the  weight  of  the  cage  is  balanced  against  the  ball/cage  forces. 
However,  large  variations  in  cage  whi.rl  velocity  are  observed  indicating 
some  element  of  instability.  The  bearing  performance  worsens  as  the  cage 
is  allowed  to  move  eccentncly  under  its  own  weight  acting  normal  to  the 
bearing  axis.  In  this  case  large  interactions  at  both  the  cylindrical  and 
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conical  parts  of  the  cage  pocket  are  observed.  This  results  in  a  large 
amount  of  noise  and  greatly  erratic  motion  of  the  cage. 

The  engine  bearing  performance  simulations  are  obtained  with  a  combined 
thrust  and  radial  load  (thrust  to  radial  load  ratio  of  4)  over  more  than 
seven  shaft  revolutions.  The  test  example  consists  of  a  100  mm  bore  ball 
bearings  operating  at  20,000  rpm  with  a  thrust  load  of  18000  N  and  a  radial 
load  of  4500  N.  The  bearing  uses  a  steel  cage  guided  on  the  inner  race. 

It  is  shown  that  for  such  an  operating  environment  the  cage  develops  appre¬ 
ciable  whirl  velocity  and  in  steady  state  a  continued  metal  contact  at  the 
guiding  land  is  established.  Hence  a  definite  possibility  of  cage  wear  at 
the  guiding  surface  is  predicted. 

Since  RAPIDREB  easily  provides  bearing  performance  simulations  over  several 
shaft  revolutions  and  the  prediction  of  steady  state  performance  is  possible 
with  great  confidence,  several  recommendations  for  future  development  may 
be  listed: 

1.  The  steady  state  performance  predictions  should  be  validated 
experimentally,  particularly  in  terms  of  the  cp?e  motion.  This 
will  greatly  strengthen  the  predictive  capability  of  RAPIDREB 
and  thereby  enhance  its  value  as  a  design  tool. 

2.  Since  solutions  over  relatively  large  time  domains  are  now 
possible  a  graphic  simulation  of  the  various  component  motions 
in  the  bearing  will  contribute  to  greater  understanding  of  the 
various  instabilities  and  the  dynamic  behavior  of  the  bearing 
as  a  whole.  This  will  essentially  require  the  development  of 
a  sophisticated  graphics  program  replacing  the  current  plot 
program. 

3.  RAPIDREB  can  be  very  useful  In  establishing  direct  correlaticns 
of  bearing  life  and  dynamic  stability  with  both  geometrical 
and  operational  variables.  Hence  it  can  be  a  very  useful  tool 
in  guiding  accelerated  life  cycle  tests  for  long  life  bearings 
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such  as  the  DMA  bearings.  A  parametric  study  to  establish  such 
correlations  and  the  development  of  guidelines  for  accelerated 
bearing  life  tests  is  therefore  recommended. 
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RDREB 
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RDREBP 
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FOR  THE  PERFORMANCE  SIMULATION  OF  RALL  BEARINGS,  THIS  PROGRAM 
IS  AN  ADVANCED  VERSION  OF  THE  BALL  REARING  TREATMENT  IN  ThE  MTI 
DYNAMICS  OF  ROLLING  ELEMENT  BEARING  (DREB)  COMPUTER  PROGRAM.  THE 
AOOEO  CAPABILITIES  *NCLUOE  — 

1.  SELECTIVE  CONSTRAINTS  OF  THE  AXIAL  AND  RADIAL  MOTION  OF 
THE  BALLS  ANO  THEREBY  ELIMINATING  HIGH  FREQUENCY  BALL/RACE 

vibrations. 

2.  DAMPING  COEFFICIENTS  AT  ALL  INTERACTIONS  WHICH  ALLOW  FOR 
EITHER  A  REALISTIC  INTERNAL  MATERIAL  DAMPING  OR  A 

fictitious  damping  to  selectively  damp  out  certain  compo¬ 
nents  of  THE  VIBRATORY  MOTION. 

3.  A  NUMBER  of  DIFFERENT  INTEGRATING  ALGORITHMS  TO  INTEGRATE 
THE  DIFFERENTIAL  EQUATIONS  OF  MOTION. 

A.  A  TREATMENT  for  a  ball  guiofd  cage. 

5.  IMPROVED  pall/racf  traction  INTEGRAL  specially  for 
solid  lurricatfo  conditions. 

6.  switch  for  simulating  a  rotating  radial  load  on  the 

BF  APING. 


RDREB  6 
■RDPEB  7 
RDREB  8 
ROREB  9 
RDREB  10 
RDREB  11 
ROREB  12 
RDREB  13 
RDREB  14 
RDREB  15 
RDREB  16 
RDREB  IT 
RDREB  18 
RDREB  19 
RDREB  20 
RDREB  21 
ROREB  22 
RDREB  23 
RDREB  24 
ROREB  25 
ROREB  26 
RDREB  27 
RDREB  28 
RDREB  29 
RDREB  30 
RDREB  31 
RDREB  32 
RDREB  33 
RDREB  34 
RDREB  35 
RDREB  36 
RDREB  37 
RDREB  38 
RDREB  39 
RDREB  40 
RDREB  41 
RDREB  A2 
RDREB  43 
RDREB  44 
RD»EB  45 
RDREB  46 
RDREB  47 
RDREB  48 
RDREB  49 
RDREB  50 
RDREB  51 
RDREB  52 
RDREB  53 
RDREB  54 
RDREB  55 
RDREB  56 
RDREB  57 
RDREB  58 
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7.  improved  computing  speeos.  proper  selection  of  cons¬ 
traints  ANO_ INTEGRATING  METHOD  CAN  EASILY  REDUCE  THE  COM¬ 
PUTATION  EFFORT  by  a  factor  of  fifteen  To  twenty  WHEN- 

COMPARED  with  the  original  orer  PROGRAM. 

The  program  -ROPER-  is  DESIGNED  FOR  ball  bearings  only  and  no 
treatment  of  roller  bearings  is  AVAILABLE. 


INPUT  DATA 


ALL  INPUT  OATA  IS  SUpplIEO  ON  DATA  CARDS  IN  COLUMNS  11  TO  70. 

The  FIRST  AND  THE  LAST  TEN  COLUMNS  ON  EACH  CARD  ARE  RESERVED  FOR 
IDENTIFICATION  purposes.  THE  USER  MAY  USE  THESE  COLUMNS  FOR  Any 
DESIRABLE  INFORMATION.  THE  INFORMATION  IS  JUST  READ  IN  AND 
PRINTED  OUT  IN  THE  LIST  OF  DATA  CARDS. 


The  PROGRAM  ALLOWS  EOR  BOTH  THE  ENGLISH  AND  SI  SYSTEMS  OF  UNITS. 
ALL  input  data  must  BE  IN  THE  SYSTEM  OF  UNIT  SELECTED  BY  THE  USER. 
The  FUNDAMENTAL  UNITS  AND  the  ABBREVIATIONS  USED  IN  THE  TWO 
SYSTEMS  ARE  DESCRIBE  IN  THE  FOLLOWING  TABLE  — 


quantity  English  unit 


SI  UNIT 


FORCE  (F) 
LENGTH  (L) 
TIME  (T) 
MASS  (MA ) 
TEMP  <TM> 


pound  force  <lbf> 
Inch  (in) 

SECOND  (S) 

pound  mass  <l«m) 
ueg  rankine  (OEG-R) 


NEWTON  (N) 

METER  (M) 

SECOND  (S) 

KILOGRAM  MASS  (KGM) 
DEG  KELVIN  (DEG-K) 


APPROPRIATE  UNITS  FOR  all  THE  PHYSICAL  QUANTITIES  CAN  BE  DERIVED 
IN  TERMS  OF  THE  AMOVE  FUNDAMENTAL  UNITS.  SIMPLY  FOR  CONVENIENCE 
all  angular  velocities  and  angular  accelerations  are  specified  in 
RPM  and  rpm/sec  RESPECTIVELY.  ALSO.  ALL  angles  are  specified  In 
oegpees. 


ALL  INPUT  VARIABLE  NAMES  BEGINNING  WITH  It  J*  K«  L*  M  AND  N  AWE 
INTEGERS  and  THE  REMAINING  are  real  NUM8ERS.  ALSOt  IT  WILL  BE 
important  that  all  Integers  are  right  adjusted  in  their  fields. 

AN  ABBREVIATION  -re-  IS  USEO  FOR  ROILING  element  in  THE  FOLLOWING 
INPUT  DESCRIPTION  and  also  IN  THE  PRINT  OUTPUT  OF  THE  PROGRAM. 

THE  VARIABLE  names*  COLUMN  fielos  and  VARIABLE  DESCRIPTIONS  are 
DOCUMENTED  BELOW  FOH  EACH  INPUT  DATA  CARD. 

CABO  1 


MODE  11-12  SWITCH  WITH  FOLLOWING  FUNCTIONS  — 

dynamic  solution  with  no 

CHANGF  IN  OPERATING  DATA  m  0 
DYNAMIC  SOLUTION  WITH 
MODIFIED  OPERATING  DATA  ON 
CAROS  13-19  ■  1 

DYNAMIC  SOLUTION  WITH  ARBI- 


ROREB  59 
RDREB  60 


RDREB  61 
RDREB  62 


RDRE 

RDRE 


63 

64 


RDREB  65 
RDREB  66 
RDREB  67 
RDREB  68 
RDREB  69 
RDREB  70 
RDREB  71 
RDREB  72 
RDREB  73 
RDREB  74 
RDREB  75 
onora  76 


ROREB  77 
RDREB  78 
RDREB  79 
RDREB  BO 
RDREB  B 1 
RDREB  B2 
RDREB  «3 
RDREB  B4 


RDPEB  «5 
RDREB  86 
RDREB  87 
RDREP  «8 
RDREB  89 
RDREB  90 
RDREB  91 
RDREB  92 
RDREB  93 
RDREB  94 
RDREB  95 
RDREB  96 
RDREB  97 
RDREB  98 
ROPER  99 
RDREB 100 
RDREBl 01 
RDREB102 
RDREB103 
RDREB 104 
RDREBl 05 
RDREB106 
RDRER107 
RDREB108 
RDREBl 09 
ROREB 1 10 
RDREB 111 
RDREB 112 
RDREBl 13 
RDREBl 14 
RDREBl 15 
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13-16 

17-19 


20-28 


29-40 


41-43 


CARO  2 


TR4RY  INITIAL  CONDITIONS  *  2 
AXIAL  AND  RADIAL  RALL  MOTION 
CONSTRAINED  BY  OUASI-STATIC 
EQUILIBRIUM  AT  EACH  STEP  ■  3 
CONSTRAINED  AXIAL  AND  RAOlAL 
BALL  motion  UNDER  PURE 
THRUST  L0A0  AND  CONSTANT 
RACE  SPEEDS  «  4 

RECORD  NUMBER  FOR  INITIAL  CONDITIONS.  IR*0  AT 
START. 

data  MONITOR  CODE.  WHEN  MODULUS  OF  (IQ-1>*  IU 
BEING  THE  STEP  NUMPEO.  ANO  MONT  VANISHES.  THE 
DIMENSIONLESS  SOLUTION  VECTOR  WILL  BE  WRITTEN 
IN  THE  master  data  FILE. 

VECTOR  OF  LENGTH  3  CONTAINING  THE  PRINT  CODES 
WITH  THE  FOLLOWING  SIGNIFICANCE  — 

NPWTU)  ■  0.  PRINT  all  ROLLING  ELEMENT  SOLUTIONS 

n,  skip  »•  polling  Element  solutions 

-N,  SKIP  all  ROLLING  ELEMENT  SOLUTIONS 

mp«t(2)«k  such  that  when  the  modulus  of  (iQ-i>» 

IQ  BEING  the  STEP  NUMBER,  and  K 
VANISHES.  THE  OUTPUT  IS  PRINTED  AT 
STEP  IQ. 

nprt < 3) *k  such  that  when  the  modulus  of  (iqp-d, 

IQP  BEING  THE  STEP  NUMBER  AT  WHICH  OUT¬ 
PUT  IS  PRINTED.  AND  K  VANISHES.  THt 

DIMENSIONLESS  solution  vector  is  also 
printed. 

each  component  HAS  a  EIELD  width  OF  3. 

VECTOR  OF  LENGTH  3  CONTAINING  THE  OPTIMIZATION 
COOES  (SEE  NOTE  5)  DEFINED  AS  — 

NOPT(l)  s  NUMBER  OF  STEPS  FOR  WHICH  INTEGRATION 

WILL  CONTINUE  WHEN  THE  TRUNCATION  tRROR 
EXCEEDS  THE  SPECIFIED  LIMIT  WITH  THE 
MINIMUM  STEP  SIZE. 

NOPT (2)  *  NUMBER  OF  STEPS  AFTER  WHICH  AN  ATTtMPT 
TO  OPTIMIZE  THE  STEP  SIZE  WILL  BE'madE. 
NOPT (3)  •  MAXIMUM  NUMBER  OF  STEPS  IN  THE  PUN. 

FACH  COMPONENT  has  A  field  WIDTH  OF  4. 

SWlTH  FOP  INTEGOATING  ALGORITHM  — 

ImET  ■  0.  FOURTH  ORDER  RUNGA-KUTTA-MERSON 

1.  FOURTH  ORDER  RUNGA-KUTTA-FEHLBERG 

2.  FOURTH  ORDER  RUNGA-KUTTA-ENGLAND 

3.  CLASSICAL  FOURTH  ORDER  RUNGA-KUTTA 

4.  FOURTH  ORDER  RALSTON-RUNGA-KUTTA 

5.  FIFTH  ORDER  KUTTA-NYSTROM 

6.  SIXTH  ORDER  MUTTA-RUNC.A-KUTTA 

7.  SECOnO  ORDER  RUNGA-KUTTA 
«.  THIRD  ORDER  RUNGA-KUTTA 

10.  A  PREDICTOR-CORRECTOR  FORMULA 

(IMET-10  RECOMMENDED  FOR  cageless 
BEARINGS  ONLY) 


SEE  NOTE  6  FOR  DISCUSSION  RELEVANT  TO  THIS  CARO. 


R0REB1I6 
RDREB1I7 
ROREB118 
RDRE8119 
RDREB120 
RDREB12I 
ROREB122 
ROREB123 
R0REB124 
R0REB125 
ROREB 1 26 
R0RE3127 
RDREB128 
RDPEP129 
ROREB130 
RDREB131 
RDREB132 
RDREB133 
ROREB 134 
RDREB135 
RDRER136 
RDREB137 
RDREB138 
ROREB 139 
RDRER140 
RDREB141 
RDREB142 
RDREB143 
PDREB144 
RDPEB145 
ROREB  U6 
RDREB147 
RDREB148 
RDREB149 
ROREB 150 
RDBEB151 
RDBEB152 
RDREB153 
ROREB 154 
ROREB 155 
ROREB 156 
RDRER157 
R0REB158 
R0REB159 
ROREB 160 
ROREB 1 61 
ROREB 162 
ROREB 1 63 
ROREB 164 

ROREB 165 
ROPER  166 
RDREH167 
ROREB 168 
RDREB169 
RDREB170 
RDREB171 
RDRE6172 
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C 


RDREB 1T3 


DINT 
0M I N 


OMAX 

TF 


m i 


35-46 

47-58 


INITIAL  DIMENSIONLESS  TIME  STEP  SIZE. 
MINIMUM  dimensionless  TIME  STEP  size. 


MAXIMUM  DIMENSIONLESS  time  step  size, 
final  DIMENSIONLESS  TIME. 


IE  IR.6T.0  AND  MOOE.EQ.O  GO  TO  CARO  21. 

IE  IR.GT.O  ANO  MOOE.EQ.l  GO  TO  CARD  13.1. 


RDREB 174 
RDRER175 
R0REB176 
RDREB177 
RDREB 178 
RDREB179 
RDREB180 
RDREB1B1 
RDREB182 


CARO  3 


KBT 

jnt 

MBR 

MHS 

ispec 


11-13  REARING  TYPE  DEFINED  as  FOLLOWS  — 
REARING  with  cage  ■  1 
CARELESS  BEARING  ■  -1 
14-16  SWITCH  FOR  TRACTION  INTEGRAL  — 


SINGLE  INTEGRAL  FOR  EHD 
CONDITIONS  *  0 

DOUBLE  INTEGRAL  FOR  SOLID 
LUR  CONDITIONS  »  1 

17-19  SWITCH  FOR  RE  AmO  RACE  MATERIALS  — 

STANDARD  OATA  *  0 

DATA  SUPPLIED  ON  CARDS 
10.1  ANO  10.2  ■  1 

20-22  SWITCH  FOR  HOUSING  AND  SHAFT  MATERIALS  — 
STANOARO  OATA  ■  0 

OATA  SUPPLIED  ON  CARDS 
10.3  AND  10.4  »  1 

23-70  PEERING  SPECIFICATION  CODE.  ANY  ALPHANUMERIC 


STRING  (A  MAXIMUM  OF  48  CHARACTERS)  BY  WHICH  ALL 
PRINT  AND  PLOT  OUTPUTS  WILL  BE  IDENTIFIED. 


CARO  4 


BOR  11-22 
OD  23-34 
SO A  35-46 

HO A  47-58 


ROHE  OR  SHAFT  OUTER  DIA  <L>. 
OUTER  RACE  OUTER  DIA  tL>. 
SHAFT  INNER  OIA  <L> . 

HOUSING  OUTER  DIA  <L>. 


CARD  5 


SF  11-34  VECTOR  OF  LENGTH  2  CONTAINING  THE  DIAMETRAL 

INTERFERENCE  <L>  AT  THE  OUTER  RACE/HOUSING  AND 
SHAFT/IN»'£R  RACE  FITS  RESPLV. 

EACH  COMPONENT  HAS  A  FIELD  WIDTH  OF  12. 


RDREB183 

RDREB 1 84 

P0REB185 

RDREB1B6 

RDREB 187 

RDREB 188 

RDREB 189 

RDREB190 

RDREB191 

RDREB 192 

RDREB 193 

RDREB 194 

R0REB195 

RDREB 196 

RDREB 197 

RDREB198 

RDREB199 

RDREB200 

RDREB201 

RDREP202 

RDREB203 

RDREB204 

RDREB205 

RDRE32P6 

RDREB207 

RDREB208 

RDREB209 

RDREB2I0 

RDREB211 

RDREB212 

RDREB213 

RDREB214 

R0REB215 

RDREB216 

RDREB21T 

RDREB218 

R0REB219 


CARD  A 

KTC 

11-16 

7RACTI0N  CODE  DEFINED  AS  FOLLOWS 

— 

SIMULATED  TRACTION  curve 

9 

0 

SHFLL  TURBO- 33  DJL 

9 

1 

5P4E  POLYPHENOL  FTHFR 

9 

2 

MIL-L-7808  OIL 

9 

3 

any  arbitrary  **odel 

9 

4 

N» 

17-22 

number  of  ROLLING  ELEMENTS. 

LU 

23-28 

SWITCH  FOR  SELECTION  DF  UNITS 

RDRER220 

RDREB221 

RDREB222 

RDREB223 

RDREB224 

R0RE«225 

RDRE«2?6 

ROREB22T 

ROREB228 

R0RE8229 


c  IN-LB  UNITS  *  1  RDREB230 
C  SI  UNITS  *  2  RDREB231 
C  NBW  29-34  NUMBER  OF  ELEMENTS  FOR  WHICH  OUTPUT  WILL  BE  RDREB232 
C  MONITORED.  UPPER  LIMIT  IS  SIX  (6).  RDREB233 
C  JBW  35-70  VECTOR  OF  LENGTH  NRW.  THE  COMPONENTS  WILL  DENOTE  R0RER234 
C  THE  INDICES  OF  THE  BEARING  ELEMENTS.  RDREB235 
C  EACH  COMPONENT  HAS  A  FIELD  WIDTH  OF  6.  RDREB236 
C  RDRE8237 
C  CARD  7  RDREB238 
C  ......  ROREB239 

C  RDBEB240 
C  OF  11-22  BALL  DIAMETER  <L>.  RDREB241 
C  DP  23-34  PITCH  DIAMETER  <L>.  RDREB242 
C  ALFA  35-46  CONTACT  ANGLE  (OEG)  OP  DIAMETRAL  PLAY  <L>.  RDREB243 

c  a  positive  value  denotes  contact  angle  and  a  rdreb244 
C  MEUATIVE  value  demotes  a  DIAMETRAL  play  equal  TO  RDRE8245 
C  THE  absolute  value  OF  ALFA.  RDREB246 
C  RCUR  < 1 )  47-58  OUTER  RACE  CUR  FACTOR.  RDPEB247 
C  RCUR (2)  59-70  INNER  RACE  CUR  FACTOR.  R0REB248 
C  RDREB249 
C  ****************•***••••**•**••••♦•••********•*•*•******************•  RDREB250 
C  *  THE  SEPIES  8  CAPOS  ARE  REQUIRED  ONLY  IN  CASE  OF  DYNAMIC  *  RDREB251 
C  *  SOLUTIONS.  MOOE.GT.O  ON  CARD  1.  *  RDREB252 
C  *•**•**♦#************#•*****•**•••••*******•**•*****•**************•*  RDREB2^3 
C  RDREB254 
C  CARD  8.1  RDRE9255 

C  - -  RDREB2S6 

C  RDREB257 
C  THIS  CARD  CONTAINS  ROLLING  ELEMENT  INERTIAL  PARAMETERS.  IF  RDREB258 
C  STANDARD  VALUES  ARE  OESIREO  THEN  SET  At L  VALUES  EQUAL  TO  ZERO.  RDRER259 
C  RDREB260 
C  XMB  11-22  PE  MASS  (MA|.  RDREB26 1 
C  XIX  23-34  RE  POLAR  MOMENT  OF  INERTIA  (MA*l**3> .  RDREB262 
C  xrr  35-46  RE  TRANSVERSE  MOMENT  OF  INERTIA  (MA*L**3>.  RDRER263 
C  RDREB264 
C  CARD  8.2  R0REB265 

C  -  RDREB266 

C  RDREB267 
C  This  CARO  CONTAINS  Cage  INERTIAL  PARAMETERS.  FOR  standard  VALUES  RDREB268 
C  OR  FOR  CAGELESS  BEARING  SET  ALL  VALUES  EQUAL  To  ZERO.  R0REB269 
C  RDREB270 
C  XmC  11-22  CAGE  MASS  <MA>.  RDREB271 
C  XKX  23-34  CAGE  POLAR  MOMENT  OF  INERTIA  (MA*L*#3>.  RDRER272 
C  XKY  35-46  CA«E  TRANSVERSE  MOMENT  OF  INERDIA  (MA*L**3) .  RDREB273 
C  RDREB274 
C  CARD  8.3  AND  8.4  RDRE8275 
C  .................  RDREP276 

C  RDREH277 
C  THE  FOLLOWING  TWO  CARDS  <I»1,2>  wIli  CONTAIN  THE  EFFECTIVE  RDREH278 

c  inertial  parameters  for  the  outer  u«n  and  inner  <i*2>  races  rdreb?t9 
C  RESPLY.  IF  THE  race  is  CONSTRAINED  r*J  ACCORDANCE  with  the  RDRER280 
C  SPECIFIED  motion  then  the  relevant  PARAMETERS  MAY  BE  SET  TO  ZERO.  RDREB2P1 
C  ALSO,  FOR  STANDARn  VALUES.  SET  A|_L  PARAMETERS  EQUAL  To  ZERO,  RDREB2°2 
C  RDREB?43 
C  XMRJII  n-22  PACE  mass  (MA).  RDRE'<2°4 
C  XJX(I)  23-34  oaCe  POLAR  MOMENT  OF  INERTIA  (MA*L**3) .  ROREB?°S 
C  X JY ( I )  35-46  PACE  TRANSVERSE  moment  OF  INERTIA  <mA*L**3>.  RDRE9286 
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RDREB287 

CARO  8.5  RDREB288 

- -  R0REB2B9 


BRDP 

11-22 

ball/race 

damping 

PATIO 

RDREB290 

BCDP 

23-34 

BALL/CAGE 

DAMPING 

RATIO 

RDPEB291 

RCOP 

35-46 

CAGE/PACE 

OAMPJnG 

PATIO 

RDREB292 

PDREB293 

R0REB294 


THE  SERIES  9  CAROS  Are  REOUIRED  ONLY  WHEN  A  CAGE  IS  PRESENT, 

(kbt.gt.o).  ir  no  cage  exists  go  to  card  io.i 


CARO  9.1 


MCR  11-16 
1HRC  17-22 
IHRC  23-28 


SWITCH  FOR  CAGE  materials  — 

STANDARD  DATA  »  0 

DATA  SUPPLIED  on  CARD  10.5  a  1 
SWITCH  FOR  RE/CaGF  CONTACT  -- 

HYDRODYNAMICS  NEGLECTED  a  0 
HYDRODYNAMICS  INCLUDED  *  1 
SWITCH  FOR  CAGE/RACE  CONTACT  — 
HYORODYNAMICS  NEGLECTED  ■  0 

hydrodynamics  included  but 

NO  SQUEEZE  FILM  DAMPING  ■  1 


*  RDRER29S 

*  RDREB296 

*  RDREB297 
RDREB298 
RDREB299 
RDREB300 
RDREB30 1 
RDREB302 
RDREB303 
RDREB304 
RDRER305 
RDREB306 
RDREB307 
RDREB308 
RDREB309 
RDREP310 


IBC  29-34 

IRC  35-40 
IRGS  41-52 


CARO  9.2 


HYDRODYNAMICS  INCLUDED 
WITH  SQUEEZE  FILM  DAMPING  ■  2 
SWITCH  FOR  RE/CAGE  CONTACT  — 

USE  HERTZIAN  SPRING  a  0 

SPRING  SUPPLIED  IN  -PSPRG-  *  1 
SWITCH  FOR  CAGE/RACE  CONTACT  — 

USE  HERTZIAN  SPRING  a  0 

SPRING  SUPPLIED  IN  -RSPRG-  a  1 
VECTOR  OF  LENGTH  2.  COMPONENTS  INDICATE  TYPE  OF 

guidance  for  The  two  lands  along  the  positive 
ano  negative  x-axis  resply  — 

OUTER  race  guidance  a  1 

INNER  RACE  GUIDANCE  a  2 

EACH  COMPONENT  has  A  EIELD  WIDTH  OE  6. 


RDREB3I1 

RDRER312 

ROREB313 

RORER314 

RDRER315 

RDRER316 

RDPER3I7 

RDREB318 

RDRER319 

RDREP320 

RDREB321 

RDREB322 

RDRER323 

RDREB324 

RDREB325 

RDREB326 

RDREB327 


COI A 

11-34 

CwC 

35-46 

CCL 

47-70 

VECTOR  OE  LENGTH  ?  CONTAINING  THE  EFFECTIVE 
OUTER  AND  INNER  DIAMETERS  CL. »  OE  CAGE. 

EACH  COMPONENT  HAS  A  FIELD  WIDTH  OE  12. 
EFFECTIVE  CAGE  WIDTH  (L>  FOR  MASS  COMPUTATION. 
VECTOR  OE  LENGTH  ?  CONTAINING  THE  EFFECTIVE 
DIAMETRAL  CLEARANCES  <L>  AT  THE  CAGE/OUTER  RACE 
ANU  CAGE/IMNER  RACF  INTERACTIONS  RESPLY  FOR  THE 
COMPUTATION  OF  CHURNING  EFFECTS. 

each  component  mas  a  field  width  of  12. 


RDREB328 
• RDREB329 
RDRER330 
RDRE8331 
ROPEP332 
RDPEB333 
RDREB334 
RDRER335 
RDREB336 
RDREB33T 


CAPOS  9.3.1  ANO  9.3.2 


RDPEB338 

R0REB339 


CARD  9.3.1  IS  ONLY  REQUIRED  IF  IPGSlD.NE.O 


RDREB340 

*  RDREB341 

*  RDREB342 

*  RDREB343 
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RDBEB344 


R0REB344 

THESE  TWO  CAROS  <r*1.2)  CONTAIN  DATA  FOR  THE  TWO  GUIDING  LANDS  RDREB345 
ALONG  THE  POSITIVF  <1*1)  AND  NEGATIVE  (1*2)  X-AXIS  RESPLY.  RDREP346 

RDREB347 

RGL(I)  11-22  GUIDING  RACE  RADIUS  <L).  RDREB348 

CRL  ( I )  23-34  GUIDING  CAGE  RAOIUS  (L).  R0REB349 

CWS(I)  35-46  LAND  WIOTH  (L)  RDREB350 

CLS(l)  47-58  DIS  OF  LAND  OUTER  EDGE  FROM  CAGE  CENTER  <L>«  R0REB351 

RDREB352 


CARD  9,4 

RCL 

11-22 

RCLH 

23-34 

PANl 

35-46 

PAN? 

47-58 

CARD  9.S 

XPOS 

11-22 

RPOS 

23-34 

APOS 

35-46 

CARO  10. 

1 

RE/CAGE  DIAMETRAL  CLEARANCE  <L> . 

RE/CAGE  MAX  LUB  FILM  (L>  .LT.  RCL  EOR  STARVED 
CONTACT.  RCLH*RCL  EOR  FULLY  flooded  CONDITIONS. 
THIS  INPUT  IS  RELEVANT  ONLY  EOR  BALL  BEARINGS. 
CONE  ANGLE  (DEG)  FOR  RE/CAGE  GUIDE  LAND. 

=  0.  FOR  NO  RE/CAGE  GUIDANCE. 

CONE  HEIGHT  <L). 

,GT.  0.  FOR  GUIDANCE  CONE  ON 
.LT.  0.  FOR  GUIOANCE  CONE  ON 


outer  diameter, 
inner  diameter. 


INITIAL  AXIAL  POSITION  ( 1_ >  OE  CAGE  MASS  CENTER 
RELATIVE  TO  THE  MEAN  AXIAL  POSITION  OF  ROLLING 
FLEMENTS. 

INITIAL  RADIAL  POSITION  <L>  OF  CAGE  MASS  CENTER 
RELATIVE  TO  ThE  CENTER  OF  THE  GUIDING  RACEWAY 
(BEARING  CENTER  IN  CASE  OF  RE  GUIDANCE) 

INITIAL  POTATION  (DEG)  OF  CAGE  AROUT  THE  X-Axjs. 


*  THIS  CARO  IS  REQUIRED  ONLY  IF  MBR«1  ON  CARD  3. 


EB 

11-22 

»***< 

RE 

ELASTIC  MODULUS 

(F/L**2>. 

POB 

23-34 

RE 

POISSON-S  RATIO. 

DEN 

35-46 

RE 

material  density 

<ma/L*#3) . 

CARD  10 

.2 

»**•< 

*  THIS  CARO  IS  REQUIRED  ONLY  IF  MBR»1  ON  CARD  3. 


ER  11-22  PACE  ELASTIC  MODULUS  (F/L**2). 
POP  23-34  RACE  POISSON-S  RATIO. 

OEmR  35-46  RACE  MATERIAL  DENSITY  (MA/L**3), 

CARO  10.3 


RDREB353 

RDRER354 

RDREB355 

RDREB356 

RDREB357 

RDREB358 

RDREB359 

RDREB360 

RDREB3M 

RDREB362 

R0RE9343 

R0REB364 

RDREB365 

RDREB366 

RDREB367 

RDREP3S8 

RDPEB349 

RDOEB370 

RDRER371 

RDREP372 

RDREB373 

RDREB374 

RDREB375 

RDRE«376 

RDREB377 

RDREB378 

RDREB379 

PDRE33B0 

RDREB3B1 

RDREB3B2 

RDREB3S3 

R0REP3B4 

RDREB385 

RDREB3B6 

RDREB3B7 

R0REB3R8 

RDREB3S9 

RDREB3R0 

RDREB391 

RDREP3R2 

RDREB393 

RDREs39h 

RDRE°395 

R0BEP3R6 

PDPES397 

RDREb39B 

RDREB399 

RDPEP400 
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*  THIS  CARO  IS  REQUIRED  ONLY  IF  MHS«1  ON  CARD  3. 


ES 
P  OS 
DENS 


11-22  SHAFT  ELASTIC  MODULUS  (F/L**2). 

23-34  SHAFT  POISSON-S  RATIO. 

35-46  SHAFT  MATERIAL  DENSITY  (MA/L#«3) . 


CARD  10.4 


*  this  card  is  required  only  if  mhs*i  on  card  3. 


eh 

pom 

DENH 


11-22  HOUSING  ELASTIC  MODULUS  (F/L**2) . 

?3-34  HOUSING  POISSON-S  RATIO. 

35-46  HOUSING  MATERIAL  DENSITY  (MA/L**3>. 


CARD  10.5 


*  THIS  CARD  IS  REQUIRED  ONLY  IF  MC»*1  ON  CARD  9.1. 


EC  1 1-22  CAUE  ELASTIC  MODULUS  <F/L**2> . 
POC  23-34  c AUE  POISSON-S  RATIO. 

DENC  35-46  CAGE  MATERIAL  DENSITY  (MA/L**3) . 

CARD  11 


JE 

11-16 

VECTOR  OF  LENGTH  ?.  COMPONENTS  DENOTE  INITIAL 
PACE  CONSTRAINTS  4i  ONG  axial  and  radial 

DIRECTIONS  (FOR  QUASI-STATIC  COMPUTATIONS)  — 
SPECIFIED  FORCE  *  0 

SPECIFIED  DISPLACEMENT  *  1 

EACH  COMPONENT  has  A  field  width  OF  3. 

JN 

17-19 

NUMBER  of  ROLLING  ELEMENTS  LAGGING  THE  ZERO  DEG 
POSITION  IN  THE  C«SE  OF  A  PARTIAL  BEARING. 

I  ROT 

20-22 

SWITCH  FOR  ROTATING  LOAD  — 

STATIONARY  RADIAL  LOAD  *  0 

ROTATING  RADIAL  LOAD  WITH 

INNER  RACE  ROTATION  *  1 

QE 

23-46 

VECTOR  OF  LENGTH  ?.  COMPONENTS  DENOTE  EITHER 

FORCE  (F)  OR  DISPLACEMENT  (L)  ALONG  THE  AXIAL  ANO 
RADUl  DIRECTIONS,  AS  PER  CONSTRAINTS  SPECIFIED 

RY  VECTOR  JF. 

faCh  COMPONENT  HAS  A  FIELD  WIDTH  OF  12. 

GAMA 

47-56 

INITIAL  MISALIGNMENT  (DEG) 

ARC 

59-70 

ANGLE  (DEG)  OVER  which  NO  RE  EXISTS  IN  ThE  CASE 

OF  PARTIAL  BEARING.  6RC=0.  FOR  FULL  BEARING. 

CARD  1? 


PPM  11-34  VECTOR  OF  LENGTH  ?.  THE  COMPONENTS  CONTAIN  T HE 
INITIAL  ANGULAR  VEIOCITIFS  (RPM)  OF  THE  OUTER 
AND  INNER  RACES  RESPLY. 


RDREB401 
RDRER402 
RDREB403 
RDRER4P4 
RDRER405 
RDRER406 
RDRER407 
RDRER408 
RDREB409 
RDREB410 
RDRE94 1 1 
RDRER412 
RDRER413 
RDRER414 
RDREB415 
RDREP416 
RDRER4 1 7 
RDRER41S 
RDRER4 1 9 
RDPER420 
RDREB421 
RDREP422 
RDRER423 
RDREP4?4 
RDREP4P5 
RDRER4?6 
RDREB427 
RDREP4P8 
RDREP4P9 
RDREP430 
RDREP431 
RDRE9432 
RDREP433 
RDRER434 
R0RER435 
RDREP436 
RDREP437 
RDRER43B 
RDREB439 
RDREB440 
RDREP44 1 
RDRER442 
RDPEP44  3 
RDPEP444 
RDREP445 
RDREB446 
RDREP447 
RDREP446 
RDREP449 
RDPER450 
RDRER451 
RDRE°452 
RDRER453 

PDREP4S4 

BDREP455 

RDREP456 

RDRER457 
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EACh  COMPONENT  was  a  FIElD  width  of  12. 
CARO  13.1  \ 


FCF1  11-58  vector  OF  length  4.  THE  COMPONENTS  DENOTE 

TRACTION  PARAMETERS  when  ktc*o.  ON  CARD  6.  OR 
WHEN  THE  actual  CONDITIONS  are  OUT  OF  LUB  MODEL 
OOUNOS.  THESE  PARAMETERS  ARE  ALSO  USED  FOR  The 

poller/flange  interactions  in  case  of  a  roller 
searing,  the  four  components  are  defined  as  — 

(1)  TPACT I  ON  COErF  AT  ZERO  SLIP. 

(2)  MAX  TRACTIOm  COEFF. 

(3!  TRACTION  COFFF  AT  infinite  slip. 

( 4 )  SLIP  (L/T*  AT  MAX  TRACTION  COEFF. 

THE  TRACTION  COEFF.  m(),  is  RELATED  TO  THE  SLIP 
VELOCITY,  V,  BY  THF  RFLATION  OF  THE  FORM  — 

MIJ  «  f  A*B*V )  *EXP  ( C •  V )  *D 

WHERE  THE  COEFFICIENTS  A,  B.  C  AND  D  ARE  DETLR- 
mINEO  BY  THE  ABOVE  FOUR  CONDITIONS.  WHEN  THE 
FOURTH  COMPONENT  OF  FCF 1  IS  SET  EOUAL  TO  ZERO, 

IT  IS  USED  AS  A  Switch  and  the  first  Three 
COMPONENTS  MUST  Thfn  CONTAIN  THE  VALUES  OF  A.  B 
and  c  DIRECTLY.  it  should  BE  noted  That  B  AND  C 
HAVE  DIMENSIONS  (T/L).  IF  THE  FOURTH  COMPONENT 
TS  NEGATIVE  THEN  ThE  VALUES  FOR  A,  B  AND  C  AWE 
STILL  CONTAINED  !►•  The  FIRST  Three  COMPONENTS, 
hO-EVER.  THE  maximum  VALUE  OF  THE  TRACTION 
COEFFICIENT  is  NO  LIMITED  to  ABS1FCF1  (4) )  , 

FACh  COMPONENT  has  A  FIELD  width  OF  12. 

CARD  13.2 


•  +<•##«#••*<•«  •««««•«•«»»#«••##«#•##«#»#«#«#•«#»••»»•  »#»»»#•»#»»»••»•»« 

*  This  Capo  IS  REQUIRED  ONLY  IF  KTc.NF.n  ON  CARD  6.  * 


IMOD 

HC 

STR 

TF« 


l 1 ”16  FILM  THICKNESS  COMPUTATION  CODE  — 

GRUBIN-S  FORMULA  »  0 

DOwSON-HIGGlNSO*'  FORMULA  «  1 

HAMROCK-DOWSON  FORMULA  a  2 

17-28  minimum  FILM  FOR  TRACTION  MODEL  BREAKDOWN  (L>. 

29-AO  STARVATION  PARAuETFR  -  DISTANCE  OF  THE  FILM 

POUNOAPY  FROM  ThE  EDGF  OF  THE  CONTACT  ZONE 
OIVIOEO  BY  THE  CONTACT  HALF  WIDTH  (SEMI-MINOR 
AXIS  IN  CASE  OF  Eli  I PT I C  CONTACT). 

M-64  VECTOR  OF  LENGTH  ?  CONTAINING  OUTER  AND  INNER 
PACE  TEMPERATURFS  (Tm>  . 
paCh  COMPONENT  mas  A  FIELD  WIDTH  OF  12. 


CARO  13.3 


*  This  CARD  IS  REQUIRED  ONLY  IF  KTc«w  ON  CARD  6. 


RDRERA58 

RDREBA59 

RDREBL60 

RDREBW61 

R0REB462 

RDREBA63 

RDREBWfeA 

RDREB465 

RDREB466 

RDREB467 

RDREP468 

RDREP469 

ROREP470 

RDREB471 

RDREP472 

HDRE0473 

PDRER474 

RQR£B475 

RDREP476 

RDREB477 

RDREP478 

RQREB479 

ROOF04PO 

rDSEbhoi 

B[)OFF4«2 

RORE04S3 

POREc4t4 

RDRFa4°5 

BDcr Q4«b 
RlyJFR4P7 
R0REP4°8 
RDREa4B9 
RDREB490 
PDPEU4R1 
BDREBW92 
RDREP4R3 
RDREB494 
R0PEH495 
R0PEH496 
RDPE=>497 
RDREH498 
PQRER499 
RDREPSOO 
RDREBSf'l 
RDPEP5,'I2 
RDREP5C3 

RDRER5f'4 
RDPEP50S 
RDREP5C'6 
RDREB5n7 
RDREP508 
RDREP5C9 
RDRER510 
RDRER51 1 
RDRERS12 
RDREB513 
RDREB51W 
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XMU  11-2?  VISCOCITY  PARAMETER  <F«T/L**2> .  R0REB515 

GAmmA  23-34  PR-VIS  PARAMETER  <l**2/F>.  R0REB516 

BETA  35-46  TE^P-VIS  PARAMETER  (1/TM).  RDREB517 

VR  47-58  CHARACTERISTIC  VELOCITY  (L/T)  IN  THE  ROLLING  R0RER518 

SPEED  EFFECT  EXPRESSION.  RDRERS19 

VN  59-70  EXPONENT  IN  THE  ROLLING  SPEEO  EFFECT  EXPRESSION.  ROREB520 

RORER521 

CARO  13.4  RDREBS22 

-  R0RER523 

RDRER524 

••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••o'***  RDRERS25 


*  THIS  CARO  IS  REQUIRED  ONLY  IF  KTC«1  or  4  ON  CARO  6.  *  RDREB526 

•••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••  RDRER527 

RDREB528 


TO 

11-2? 

INLET  TEMPERATURE  (TM). 

R0REB529 

XmuO 

23-34 

INLET  VISCOCITY  (F*T/L*«2). 

RORERS30 

GAM 

35-46 

PR-VIS  COEFFICIENT  (l**2/F) . 

RDPEB531 

PET 

47-58 

TEMp-viS  COEFFICIENT  (TM). 

RDREB532 

FK 

59-70 

THtRMAL  CONDUCTIVITY  (F/T/TM). 

RDRE8533 

CARD  1 4 

1CHR 

11-16 

SWITCH  FOR  CHURNING  FFFECTS  — 

RDRER534 

RORERS35 

RDREB536 

RDREB537 

CRO 

17-28 

CHURNING  NEGLECTED  * 

CHURNING  CONSIDFREO  « 

EFFECTIVE  VISCOCITY  is 

DERIVED  FROM  LtJa  -OPEL 
(KTC.GT.O  ON  Capo  6)  * 

EFFECTIVE  LUB  DENSITY  (MA/L**3). 

0 

1 

2 

RDRER53B 

RDREPS39 

RDRER540 

RDREP541 

RDRER542 

RDRER543 

vise 

29-40 

EFFECTIVE  VISCOCITy  <E*T/L**2). 

IF  I CHR*2  ON 

RDREP544 

CARD  15. 

1 

THIS  CARD  THEN  EPF  VIS  ■  VISC*(LUB  MODEL  VIS). 

RDRER545 
RDRER546 
RDRER547 
RDRE«548 
**  RDRE°549 

*  This  CARD  IS  REQUIRED  only  if  a  cagf  is  present. 

*  RDREB550 

*  krt.gt.o 

ON  CARD 

3. 

*  ROREP551 

••  RDPER552 

RDRER553 

FCF 2  11-58  VECTOR  OF  LENTTH  4  CONTAINING  THE  TRACTION  PARA-  RDRER554 

METERS  FOR  RE/CAGF  INTERACTION.  THE  PARAMETERS  RD«ER555 
HAVE  SAME  DEFINITIONS  AS  FOR  FCF1  ON  CARO  13.1,  R0RER556 
FACH  COMPONENT  HAS  A  FIELD  WIDTH  OF  12.  RDREB557 

HCRP  59-70  CRITICAL  FILM  (L)  FOR  RE/CAGE  CONTACT.  RORER558 

RDRER559 

CARO  15.2  RDREB560 

- -  RDREB561 

RDREB562 

••«••»•••*••••»•••••••»••••••••••••••••••••••••••••••••••••••••••••••  RDRER563 

*  THIS  CARD  IS  REQUIRED  ONLY  IF  A  C AGF  IS  GUIDED  ON  THE  RACE.  *  R0RER564 

*  KRT.GT.O  ON  CARD  3  AND  IRGS(i)*IRGS(2)  .GT.O  ON  CARO  9.1.  •  RDRER545 

***•***•«••••**•****«****•*••••*•*••*•***•*•***•*•*••****♦*******•***  R0RE85A6 

RDREBSA7 

FCF3  11-58  vECrno  OF  LENGTH  4  CONTAINING  THE  TRACTION  PARA-  R0®EP568 

mETERS  OF  RACE/CAOr  INTERACTION.  THE  PARAMETERS  RDRER569 

have  Same  DEFINITIONS  AS  FOR  FCF 1  ON  CARD  13.1.  RDREB570 
FACh  COMPONENT  HAS  A  FIELD  WIDTH  OF  12.  RDRER571 


-100- 


ooooo  ooooooooooooooooo  oo  r>  o  n  o  r>  nooooooooo  on  oooooooooooonooo 


HCC®  59-70  CRITICAL  FILM  <L>  FOR  RACE/CAGE  CONTACT. 
CARO  16 


*  this  caro  is  required  only  if  a  cagf  is  present.  * 

*  kbt.gt.o  on  capo  i.  * 


CFOR 


RTC 

THC 


11-46  VECTOR  OF  LENGTH  3  CONTAINING  THE  x«  y  AND  z 
COMPONENTS  (INERTIAL  FRAME)  OF  ANY  EXTERNAL 
FORCE  (F)  vector  acting  on  the  cage*  such  as 

LUBRICANT  JET  FORCF. 

EACH  COMPONENT  has  a  field  width  OF  12. 

47-58  DAOIAL  COORDINATE  (L)  RELATIVE  TO  THE  BEARING 
CENTER  AT  WHICH  The  ABOVE  FORCE  IS  APPLIED. 
59-70  ORBITAL  COORDINATE  (DEG)  DEFINING  THE  POSITION 
VECTOR  LOCATING  Thf  POINT  OF  APPLICATION  OF 

cfor  in  the  inertial  frame. 


CARD  17 


*  this  card  is  required  only  if  ktc»o  on  card  6.  kbt.gt.o  on  * 

•  CARD  3  ANO  I H9C*  TmRC  .GT .  0  ON  CAR'-'  B.l.  * 


VISP  11-2?  FFF  LUB  VIS  <F*T/l**?>  FOR  RE/CAGE  INTERACTION. 

a  zero  may  be  specified  if  ihbcso  on  card  9.1. 

VI SR  23-34  FFF  hjb  VIS  <E*T/l**2>  FOR  RACE/CAGE  INTERACTION. 

a  zero  may  be  specified  if  ihrc«o  on  card  9.1. 

CARD  18 


«  this  capo  is  requireo  only  if  mode.ge.o  on  card  1. 


the  following  card  specifies  race  constraints  for  dynamic 

SOLUTIONS.  IRF  AND  IRM  ARE  TWO  0I“FNST ONAL  ARRAYS  (2»3) .  THt 
FIRST  INDEX  DENOTES  THE  RACE  <OUTER=l»  INNER«2)  anD  The  SECOND 
INDEX  DENOTES  THE  Three  COMPONENTS  OF  THE  INERTIAL  TRIAD  — 


IRF  11-28  LINEAR  CONSTRAINT  VECTOR  OF  LENGTH  6.  THE 
COMPONENTS  ARE  DEFINED  AS  — 

SPFCIFIEO  FORCE  a  0 

SPECIFIED  DISPLACEMENT  a  1 

taCh  COMPONENT  HAS  A  FIELD  WIDTH  OF  3. 

I R’l  29-46  ANGULAR  CONSTRAINT  VECTOR  OF  LENGTH  6.  THE 
COMPONENTS  ARE  DEFINED  AS  — 

SPECIFIED  MO»E‘<T  a  0 

SPECIFIED  ANG  DISPLACEMENT  =  1 

FACh  COMPONENT  has  A  FIELD  WIDTH  OF  3. 


RDRER572 

RDREB573 

RDREB574 

RDREB575 

RDRE8576 

RDREP577 

RDREB578 

RDREP579 

RDRE5580 

RDREPSei 

RDRE8582 

RDREB5B3 

R0REBS84 

RDREB5«5 
RDREPS86 
RDREB5P7 
RD»E85«S 
RDRERSP9 
RDPE°S90 
R0REB59 1 
R0REP592 
RDREP593 
RDRE9594 
PDRE8595 
RDREP596 
RDREP597 
RDRES598 
RDREB599 
RDRE8600 
RDREPSOl 

rdrepg^z 

RDREP6P3 

RDPEP6P4 

RDREP605 

RDREP606 
RDRER6C  7 
RDREP6D8 
RDREPt 09 
RDRE8610 
RDREP61 1 
RDRE8612 
RDPEB613 
RDREP614 
RDREP615 
RDREP616 
RDREP617 
RDREP618 
RDREP61 9 
RDPE8620 
RDPEB621 
RDREP622 
PQREP623 
PDREn6?4 
PDCEB6?5 
PDREB626 
RDRE8627 
PDPrP620 
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CARO  19 


this  caro  is  required  only  if  mode.ge.o  on  card  i. 


6  11-46  VECTOR  OF  LENGTH  3  CONTAINING  THE  COMPONENTS  OF 

THE  ACCELERATION  <L/T*«2)  DUE  TO  GRAVITY  VECTOR 
IN  THE  INERTIAL  FRAME.  EACH  COMPONENT  HAS  A 
FIELO  WIDTH  OF  12. 


THE  FOLLOWING  SERIES  20  CARDS  ARE  REQUIRED  ONLY  IF  M0DE®2  ON 
CARD  1. 


CARD  20.1 


Tl  11-25  INITIAL  DIMENSIONLESS  TIME. 
CARDS  20.2  TO  CARD  20.N 


X  11-70  THE  DIMENSIONLESS  SOLUTION  VECTOR  OF  LENGTH 

(30*12*NB»6*FBT/A8S(<RTn  WHICH  WILL  8E  USED  AS 
INITIAL  CONDITIONS.  EACH  CARD  WILL  HAVE  FOUR 
COMPONENTS.  EACH  one  HAVING  A  FIELD  WIDTH  OF  15. 

CARDS  21  TO  . 


ANY  INPUT  DATA  REQUIRED  BY  THE  USER  ACCESSIBLE  ROUTINES  -ALOAD-, 
-PSPRG-  AMD  -RSPRO-  SHOULD  BE  SUPPLIED  IN  ORDER.  (SEE  DISCUSSION 
OF  SPECIAL  ROUT  I  NFS  BELOW). 


OUTPUT  DATA 
ssttott.a: 

OUTPUT  SOLUTIONS  may  be  EITHER  PRINTED,  PLOTTED  OR  BOTH.  AS  DIS¬ 
CUSSED  UNDER  THE  INPUT  DESCRIPTIONS  SEVERAL  DIFFERENT  OPTIONS  FOR 
MONITORING  THE  OUTPUT  APE  AVAILAPLF.  ALL  THE  OUTPUT  IS  VERY  WELL 
DOCUMENTED  IN  ORDER  TO  MAKE  IT  SELF-EXPLANATORY.  ALL  DIMENSIONAL 
VARIABLES  ARE  PRInTEo  WITH  PROPER  UNITS. 


file  systfm 


THE  PROGRAM  HAS  FULL  9FSTART  CAPABILITIES  AND  IT  CAN  MONITOR  IhE 
MOTION  OF  A  TOTAL  OF  six  BEARING  ElFmEnTS  INCLUDING  THE  CAGE  *ND 

the  raceways,  the  Element  selection  is  dfscribed  under  the  input 

DESCRIRTION  AND  SOMt.  DETAILS  OF  TrtF  FILING  SYSTEM  ARE  DISCUSSLO 
BELOW. 

FILE  esol 


RDRE8629 
RDREB630 
R0REB63 1 

*  RDREB632 

*  ROREB633 

«  RDREB634 

RDREB635 

RDREB636 

RDREB637 

R0REB638 

RDREBS39 

RDREB640 

RDRER641 

*  RDPER642 

*  RDREB643 

*  RDRER644 

*  RDPEB845 
RDREB646 
PDREB647 
R0RE«648 
RDRER649 
ROREB6SO 
RDREB651 
RDREB652 
RDRER653 
RDREB654 
RDREB6S5 
RDREB656 
RDREB6S7 
BDREB658 
RDRER659 
RDREB660 
RDREB661 
RDRE«662 
RDREB663 
RDREB664 
RDREB665 
RDREB566 
R0REB667 
RDREP688 
RD»EB6G9 
RDPEB670 
RDREB671 
RDREB672 
RDREB673 
RDREB674 
RDREB675 
ROBERT 
RDREPS77 
RDPEB878 
RDBER6T9 
RDREBftRO 
RDREB8P1 
PDBEB892 
RDREB683 
R0REB6«4 
RDREB685 
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C 


R0REB6R6 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  IS  the  MASTER  Data  FILE  WHICH  CONTAINS  ALL  THE  BEARING 
INFORMATION  AND  SOLUTION  VECTORS  at  SELECTED  time  steps.  This 
file  MUST  re  SAVEn  If  the  execution  has  to  be  restarted  LATER. 


file  FSOL 


AT  THE  END  OF  A  SUCCESSFUL  RUN  THE  LAST  SOLUTION  VECTOR  IS 
WRITTEN  IN  this  file  IN  A  FORMAT  SlJCH  THAT  THE  DATA  COULD  BE 

punched  on  cards  anu  directly  used  under  the  arbitrary  initial 

CONDITIONS  OPTION.  THIS  FILE  IS  COMPLETELY  REWRITTEN  AT  THE 
END  OF  EACH  EXECUTION  OF  THE  PROGRAM. 


FILES  GSOL1  TO  GSOLP 


these  six  files  store  all  solutions  for  the  six  preselected  bear¬ 
ing  elements.  if  the  number  OF  ELEMENTS  SELECTED  IS  LESS  THAN  SIX 
then  the  first  file  numbers  are  filiEd.  in  other  words  if  there 

are  ONLY  three  ELEMENTS  MONITORED  THEN  THE  SOLUTIONS  ARE  CONTAIN¬ 
ED  IN  FILES  GSOL1  TO  C.SOL3  AND  ThE  FILES  GSOL4  TO  GS0L6  REMAIN 

Empty,  when  a  complete  case  consists  of  several  executions  of 

THE  PROGRAM,  Then  IT  WILL  BE  NECESSARY  TO  ATTACH  ALL  THE  FILES 
USEO  WITH  THE  PROPER  WRITE  PERMISSION  SO  THAT  SOLUTIONS  FOR  TH£ 
ENTIRE  TIMF  range  may  BE  STORED  FOR  PLOTTING  PURPOSES. 

FILE  HSOL 


The  over  all  BEARTNR  PARAMETERS,  which  consists  OF  L0AD*SLIP« 

TRACT  I ON*SL IP.  ANO  TORQUE  VARIATIONS.  ARE  STORED  IN  THIS  FILE  AS 
A  FUNCTION  OF  TIMF, 

ALL  DATA  FILES  CONTAIN  APPROPRIATE  IDENTIFICATION  CODES,  WHICH  ARE 
CHECKEO  whenever  THt  FILES  ARE  ATTACHED.  THE  FILES  PRODUCED  »ITH 
TWO  DIFFERENT  RUNS.  THEREFORE.  MAY  NOT  BE  INTERCHANGED. 


SPECIAL  ROUTINES 

H>t*:nsK«K» 

the  following  special  routines  are  available  to  the  user  in  okder 
to  CLOSELY  simulate  THE  APPLICATION  UNDER  QUESTION  — 

1.  APLOAP  -  APPLIED  LOAD 

2.  PSPRG  -  ARBITRARY  SPRING  AT  RAlL/CAGE  CONTACT 

3.  PSPRG  -  ARBITRARY  SPRING  AT  RACE/CAGE  CONTACT 

SEE  DOCUMENTATION  SUPPLIED  WITH  THE  ABOVE  ROUTINES  FOR  NOTES  On 
programming  The  FEATURES  OF  A  GIVEN  APPLICATION  IN  the  ABOVE 
SUBROUTINES. 


EXTERNALS  CALLED 

■  imictistiiiii 

THE  FOLLOWING  LIST  OF  SUBROUTINES  INCLUDES  THE  ROUTINES  CALLED 

directly  by  -BnRFo-  and  also  those  callfd  by  the  routines 
initially  called  by  -rdpep-  — 


RDREB6B7 
RDREB6B8 
RDREB6P9 
RD»EBF>90 
R0REP691 
RDREB692 
RDREB693 
RDRER694 
RDREB695 
RDPEB896 
PDREB697 
R0REB698 
RDREB699 
RDRER700 
ROPEBTOI 
PDRER702 
RDREBT03 
RDREB704 
RDPEB705 
RDREB706 
PDREB7A7 
RDREB708 
RDRER709 
RDRER7] 0 
PDREB7U 
RDPEB712 
RDREP713 
R0»EB714 

RDRE»7)5 
RDREP716 
RDREB717 
RDREB718 
RDREB7’ 9 
RDREP720 
RDREP721 
RDREB722 
RDREB723 
RDREB724 
RDPEB725 
RDREB7?6 
RDRER727 
RDREP728 
R0REB7?9 
RDREP730 
RDREB731 
RDREB732 
RDREP733 
R0REP734 
RDREB735 
RDREP736 
BQRE«7?7 
RDREb7?8 
P0OEP719 
RDSEP740 
RDRE°7<.1 
R0REP742 
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c 

RDREB743 

c 

1. 

ABM 

m 

PREDICTOR-CORRECTOR  ALGORITHM 

RDREB744 

c 

2. 

ACCN 

- 

ANGULAR  ACCELERATION  transformation 

RDREB745 

c 

3. 

ARTAN 

m 

MOOIFIED  arctan 

RDREB746 

c 

4. 

PCCON 

- 

pall/cage  interaction 

R0REB747 

c 

5. 

9GEOM 

■k 

bearing  geometry 

RDREB748 

c 

6. 

BLOAO 

•k 

rall/race  loao 

RDREB749 

c 

7. 

BT9AC 

a* 

rall/pace  traction 

RDREB750 

c 

8. 

CELIL 

m 

ELLIPTIC  INTEGRALS 

RDREB751 

c 

9. 

CHURN 

m 

CHURN iHG  MOMENTS 

RDREB752 

c 

10. 

crcon 

m 

CAGE/RACF  INTERACTION 

RDREB753 

c 

11. 

CUPP 

a* 

curvature  function  for  elastic  contact 

RDREB754 

c 

12. 

DERIV 

- 

COMPUTATION  OF  DERIVATIVES 

RDRER755 

c 

13. 

ELCON 

m 

elasticity  solution 

R0REB756 

c 

14. 

FCOEF 

traction  coefficient 

R0REB757 

c 

15. 

FCT 

- 

CURVATURE  function  for  elastic  contact 

RDREB758 

c 

16. 

FCTT 

- 

cage/pace  contact  coordinate  function 

PDREB759 

c 

17. 

FLB 

• 

long  reaping  FUNCTION 

RDREP760 

c 

18. 

FLIFE 

- 

FATIGUE  LIFE  COMPUTATION 

RDR£B76 1 

c 

19. 

FPHI 

«• 

pale  AXIMUTH 

RDREB762 

c 

20. 

FR1C 

•» 

friction  FORCE  FUNCTION 

RDRER7A3 

c 

21. 

FSB 

«• 

short  BEARING  function 

RDREB764 

c 

22. 

FTPS 

m 

TRACTION  PARAMETER  FUNCTION 

RDREP765 

c 

23. 

FIJNC 

• 

INTEGRAND  for  TRACTION  integral 

RDREB766 

c 

25. 

hydro 

hydrodynamics  fop  short  and  long  bearing  models 

RDREB767 

c 

26. 

I  ANG 

• 

INITIAL  ANGULAR  VELOCITIES 

RDREB7A8 

c 

27. 

PCOFFF 

polynomial  COEFFICIENTS  for  PREDICTOR-CORRECTOR 

RDRER769 

c 

28. 

PHITS 

m 

thermal  and  SIDE  leakage  FACTORS 

RDPEP770 

c 

29, 

PLOAD 

• 

PALmgpEN  FORMULA 

RDREB771 

c 

30. 

POUT 

• 

PRINT  OUTPUT 

RDREB772 

c 

31. 

prvit 

• 

PRESSURE- V I SCOC IT Y-TEmPERATURE  RELATION 

RDREBT73 

c 

32. 

OGIO 

m 

TE <Th  ORDER  GAUSSIAN  QUADRATURE  INTEGRATION 

RDREB774 

c 

33. 

ostat 

m 

DUASI-STATIC  equilibrium 

RDREB775 

c 

34. 

9KM4 

m 

explicit  INTEGRATING  ALGORITHMS 

RDREP776 

c 

35. 

rollf 

m 

POLLING  FUNCTION 

RDBEB777 

c 

36. 

RTMT 

- 

POOT  EVALUATION  FOP  A  GIVEN  FUNCTION 

RDREB778 

c 

37. 

SUP 

- 

SLIP  VELOCITY  BETWEEN  TWO  INTERACTING  BODIES 

RDREP779 

c 

36. 

STEOM 

- 

static  equilibrium 

RDREB760 

c 

39. 

TPAC 

- 

TRACTION  ON  INCREMENTAL  AREA 

RDPE87P 1 

c 

40. 

TRCF 

• 

TRACTION  DATA  FOR  MINERAL  OILS 

RDREB782 

c 

41. 

TRFM 

- 

TRANSFORMATION  matrix 

RDREB783 

c 

42. 

TRMO0 

- 

FHO  TRACTION  MODELS 

RDPEB7R4 

c 

43. 

TRSLP 

- 

COEFFICIENTS  FOR  STMPlIFIEO  TRACTION  RELATION 

R0REB785 

c 

44. 

VEUC 

- 

angular  velocity  TRANSFORMATION 

RDREB786 

c 

45. 

VPDT 

- 

VECTOR  PRODUCT 

RDRERTR7 

c 

46. 

7  TNT 

- 

LINEAR  INTERPOLATION 

RDREB788 

c 

47. 

AUOAO 

- 

APPLIED  load 

RDREPTP9 

c 

48. 

PSPRG 

- 

arbitrary  spring  at  PALL/CAGE  contact 

RDREB790 

c 

49. 

RSPRG 

• 

ARBITRARY  SPRING  AT vHACE/CAGE  contact 

RDREP791 

c 

RDREP792 

c 

RDREPT93 

c 

notes 

RDREPT94 

c 

mu 

RDREP795 

c 

RDRER796 

c 

<1> 

A  maximum  of  “0  ROLLING  ELFme»TS  are  allowed.  if  The  actual 

RDPE°T97 

c 

MUMRfO 

EXCEEDS  40*  THEN  ALL  n I MFNS I ON  DECLARATIONS*  INCLUD- 

RDREB798 

c 

ING  THOSE 

IN  Thf  COMMON  blocks  •‘ILL  HAVE  TO  BE  CMANGEO 

RDRER799 

c 

ACCORDINGLY. 

RDRERPOO 

c 

RDREB801 
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(2)  THE  PROGRAM  WITH  ALL  THE  ROUTINES  USES  A  CENTRAL  MEMORT  OF 
ABOUT  1 580008  WORDS, 

(3)  THIS  PROGRAM  has  SEEN  DEVELOPED  ON  THE  CDC  SYSTEM  WHICH 
EMPLOYS  60  PIT  WORDS.  ON  OTHfR  SYSTEMS  WHERE  SHORTER  WORDS 
ARE  USED  IT  MAy  RE  NECESSARY  TO  USE  DOUBLE  PRECISION. 

(4)  FOP  THE  PURPOSE  OF  ROOT  FINDING  AND  TEST  AGAINST  ZEROS  The 

program  uses  tolerances  which  reside  in  the  array  tolr  In 
the  common  pluck  tcoma.  the  aroay  has  three  components. 

THE  FIRST  AND  SECOND  COMPONENTS  ARE  TOLERANCES  FOR  ABSOLUTE 
AND  RELATIVE  ZEROS  RESPLY.  The  THIRO  COMPONENT  IS  A 
TOLERANCE  ON  ACCELERATION  COMPONENTS.  if  any  ext  force 
has  AN  ABSOLUTE  value  LESS  Than  THIS  TOLERANCE  THEN  IT  IS 
set  equal  to  zero,  proper  values  to  the  array  tolr  are 

ASSIGNED  In  the  PEGINNING  OF  THE  MAIN  PROGRAM  -DREB-. 

(5)  TIME  STEP  SIZE  OPTIMIZATION  is  PROVIDEO  BY  NOPT  VECTOR  and 
THE  VALUE  OF  Tol  in  -RDREB-.  TOL  SHOULD  BE  1  TO  2  ORDERS 

of  magnituof  greater  than  thf  tolerance  on  relative  zeros 

DISCUSSED  APOVE.  IF  THE  TRUNCATION  AT  ANY  STEP  EXCEEDS  TOL 
THEN  THE  STEP  SIZE  IS  SUCCESSIVELY  HALVED  UNTIL  EITHER  The 
TOLERANCE  CRETERION  IS  MET  OR  thf  MIN  STEP  SIZE.  SPECIFIED 
ON  CARD  2.  IS  REACHED.  IF  Thf  TOLERANCE  CRETERION  IS  NOT 
MET  FOR  NOPT(l)  CONSECUTIVE  STEPS  THEN  The  EXECUTION  IS  TER¬ 
MINATED.  A  TYPICAL  VALUE  FOR  NOPTtl)  MAy  BE  ABOUT  3  TO  6. 

IF  The  TOLERANCE  CRITERION  IS  mFT  FOR  NOPT<?>  SUCCESSIVE 
STEPS  Them  am  ATTEMPT  IS  made  To  SUCCESSIVELY  DOUBLE  THh 
RTEP  SIZE  UmTIl  EITHER  THE  “A*  STEP  SIZE  OR  am  UNACCEPTABLE 
TRUNCATION  RRHOR  IS  REACHED.  T-tF  OPTIMUM  STEP  IS  THEN 
SELECTED  FOR  THE  NEXT  NOPT (2)  STFPS.  A  TYPICAL  VALUE  EUR 
MOPT(2)  may  BE  10.  IF  NO  OPTIMIZATION  IS  REQUIRED  THEN 
»'OPT(2)  MAY  BE  SET  TO  A  VALUE  GREaTFR  than  the  MAX  number 
OF  STEPS  In  the  RIJN,  WHICH  IS  SPFCIFIED  AS  nOPT<3>.  THt 
maximum  mumper  OF  STEPS  SHOIIL"  °f  selected  on  The  BASIS  OF 
the  expected  computing  time  per  step  and  the  permissible 

T I MF  LIMIT  ON  THE  JOB.  ADEQUATE  ALLOWANCE  MUST  BE  MADE  IN 
ORDER  TO  PRFVENT  -TIME  LIMIT-  ERROR*  WHICH  ON  SOME  SYSTEMS 
may  RESULT  IN  IMPROPER  CLOSUpf  or  DATA  FILES  AND  HENCE  A 
SUBSTANTIAL  PURTION  OF  the  output  MAY  be  LOST. 

(6)  FOR  MOST  BEARINGS  A  TIME  STEP  SIZE  OF  ABOUT  10  MICROSECS 
WILL  BE  ADEQUATE  FOB  CONVERGENCE  TO  WITHIN  REASONABLE 
TOLERANCE  LIMITS.  A  TRTAL  RUN  WILL  PRINT  THE  RELEVANT  Time 
SCALE  WHICH  CAN  BE  USED  TO  ESTIMATE  THE  DIMENSIONLESS  STEP 
SIZE.  THE  “IN  and  max  STEP  StZES.  SPECIFIED  AS  DMIN  AND 
Dma x  RESPLY  ON  CARO  2*  MAY  BE  ESTIMATED  BY  THE  TRANSIENT 
DETAILS  REQUIRED  IN  THE  SIMULATIONS.  THE  INITIAL  STEP  ' 
SIZE.  DINT,  MUST  BE  BETWEE*1  D*'IN  AND  OMAx.  IF  DINT*0mIn« 

Dm A X  Then  no  STFP  OPTIMIZATION  w ILL  BE  PERFORMED  AND  UNDER 
SUCH  A  CONDITION  NOPT ( 2 )  O'.  CARD  \  SHOULD  BE  SET  EQUAL  TO 

A  VALUE  GREATER  Than  THE  max  NUMPfR  OF  STEPS  IN  THE  RUN l 


PBADFEp  K.  GUPTA 

VERSION  PAPIDUEP.O  J"LK  1981 


RDREBB02 

RDREB803 

RORER804 

RDREB805 
ROREB806 
RDREP807 
RDREB808 
RDRE8809 
PDREPBl 0 
R0REP81 1 
RDPEB812 
RDREP813 
R0RER814 
PDREB815 
RDREP816 
RDREB817 
RDREB818 
RDREBB19 
RDREB820 
RDREB821 
RDPEBR22 
RDREB8P3 
RDREP824 
R0REP825 
RDREB826 
RDREB827 
RDREP828 

RDRE&830 
R0°£R8  3 1 
RDREQP32 
RDRE°833 

RDBEPB34 

RDREB835 

RDREB836 

RDREB937 

RDREB»38 

RDREP839 

RDREB840 

RDREB8-1 

RDREP842 

RDREB843 

RDREB844 

RDRER845 

RDREP846 

RDREB847 

R0REP84B 

RDRE«B49 

RDREB«50 

RDRE«851 

PDRER852 

RDREQBB3 

RQREB854 

RDREP855 

RDREh856 

ROPE(1B57 

ropebb=8 

•RDREurB9 

RDREB880 
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PROGRAM  rdrebp 


rapid  simulation  of  the 


DYNAMICS  OF  ROLLING  ELEMENT  REAPINGS  (PLOT  PROGRAM) 

sssssaas  »*  ussMi  ■■»■■■  siniiti  >iix«atiiBitn 

-PLOT  OUTPUT  For  THE  real  TI'af  PERFORMANCE  SIMULATION- 

«Y 


PRAOEEP  K.  GUPTA 

MECHANICAL  TECHNOLOGY  incorporated 
068  AL8ANY-SMAKER  ROAD 
LATHAM*  NF.^  YORK 


RDRERP  6 
RDRERP  7 
RDRERP  8 
RDRERP  9 
ROREBPIO 
RDRERP II 
R0RERP12 
RDRERP 13 
RDRERP 14 
RDRERP15 
RDRE9P16 
RDRERP17 
RDRERP18 
RDRERP19 
RDREPP20 
RDREBP21 
RDRERP22 
RDRERP23 
RDRERP24 
R0RERP25 


U.  S.  A.  RDRERP26 

RDRERP27 

SPONSORED  PY  R0RERP28 

wrighT -patterson  air  force  rase  rd»erp29 

OHIO  RDRERP30 

CONTHACT  NO.  F33615-B0-C-5152  RDPERP31 

RDREPP32 

R0REBP33 


ALL  OUTPUT  SOLUTIONS  GENERATED  By  THE  DYNAMICS  OF  ROLLING  ELEMENT  RDRERP34 
all  OUTPUT  SOLUTION  GENERATED  BY  THE  PROGRAM  -RDREB-  (RAPID  SlMU-  RDREBP35 
LATION  OF  THE  DYNAMICS  OF  ROLLING  ELEMENT  BEARINGS)  ARE  PLOTTED  BY  RDREBP36 
THIS  PROGRAM.  ThF  Data  FILE  CONTAINING  THE  SOLUTIONS  FOR  A  GIVEN  RDREBP37 
BEARING  ELEMENT  IS  attaCHEO  BEFORE  EXECUTING  -RDREBP-.  TO  MINI-  R0REBP38 
MIZE  The  COPE  REQUIREMENT  this  MaSTFR  DATA  FILE  IS  SUBDIVIDED  Into  RDREBP39 
SEVERAL  local  FILES  SUCH  THAT  EACH  LOCAL  FILE  CONTAINS  ALL  THE  RDREBPAO 
DATA  FOR  A  PARTICULAR  PLOT.  THE  LOCAL  FILES  ARE  THEN  PLOTTED  IN  RDREBP41 
SEQUENCE.  AS  DISCUSSED  BELOW  THE  PROGRAM  HAS  FULL  CAPABILITIES  OF  RDREBP42 
AUTOMATIC  SCALING  TO  any  DESIRED  DEGREE.  RQREBP43 

RDREBP44 


INPUT  description 


ALL  DATA  IS  SUPPLIED  IN  COLUMNS  11  To  TO.  THE  FIRST  TEN  AND  The 
last  ten  COLUMNS  are  AVAILABLE  fop  ANY  IDENTIFICATION  purpose; 
The  FOLLOWING  DESCRIPTION  SUMMARIZES  The  input  variable  names* 
THEIR  FIELDS  ON  the  DATA  CARDS  AnO  ThF  RELEVANT  DESCRIPTION*  IN 
ORDER.  VARIABLES  BCGINMING  WITH  Thp  (.ETTFRS  I*  J.  K*  L»  M  OR  N 
ARE  INTEGERS  And  all  OTHERS  ARE  REAL  NUMBERS.  THF  FORMAT  OF  All 
RFAl  NijmBFRS  IS  E12*5  AND  THAT  FDR  THE  INTEGERS  IS  DETERMINED  BY 
The  SPECIEICO  FIElO  LENGTH  ANO  Th£  nu-qeR  IS  RIGHT  JUSTIFIED  IN 
ITS  FIELO. 

CARD  1 


RDRERP45 

RDREBPA6 

R0REBP47 

RDREBP48 

RDREBP49 

RDREBP50 

R0REBP51 

R0RERP52 

RDREBP53 

R0REBpe>4 

RDREBPS5 

R0REBPS6 

R0REBP57 

R0REBP58 


ooonooooooooooooooooo onoooo  no  noooononooooooooooooooooo  oo 


Ml™ 


c 


I  TP 
JTP 


JR  I 
JRF 
I  PUT 


11-12  SWITCH  FOR  PLOTTING  TITLF  PAGE  ~ 

I Tp  *0  puot  title  page. 

*1  suppress  title  page  plotting. 

13-16  NUMBER  OF  POINTS  TO  PE  SKIPPED  IN  THE  PLOT 

JTPsO  WILL  PLOT  Au  POINTS.  THE  TOTAL  NUMBER  OF 
POINTS  TO  BE  PLOTTED  MUST  BE  .LE.  3000 
17-23  RECORD  NUMBER  FROM  WHICH  PLOT  SHOULD  START. 

24-30  RECORD  NUMBER  AT  WHICH  PLOT  SHOULD  END. 

31-70  A  VECTOR  WHOSE  LENGTH  IS  EQUAL  TO  THE  NUMBER  OF 
PLOTS  FOR  THE  BEARING  ELEMENT  UNOER  CONSIDERAT¬ 
ION.  there  could  pe  a  maximum  of  ten  plots, 
the  VALUE  ASSIGNED  TO  each  COMPONENT  has  The  FOL¬ 
LOWING  SIGNIFICancf  — 
i put < i )  «o  itm  plot  is  suppressed. 

si  AUTOMATIC  SCALING  IS  USED  ON  THt  ITH 
PLOT. 

»2  SCALE  INFORMATION  IS  SUPPLIED  IN  THE 
FOLLOWING  CARDS. 

EACH  COMPONENT  has  A  FIELD  WIDTH  OF  4. 


IF  NONE  OF  THE  COMPONENTS  OF  IPLT  APE  SET  EQUAL  TO  2  THEN  NO  AD¬ 
DITIONAL  DATA  CAPOS  ARE  NECESSARY.  IF  THE  VALUE  2  IS  ASSIGNED  TO 
ANY  COMPONENT  THEN  CaPOS  2  TO  5  MUST  PE  SUPPLIED  FOR  EACH  COMPON¬ 
ENT  OF  IPLT  (*2>  IN  ORDER. 


CARD  2 


ISCX(I) 


XIX(I) 
XFX ( I  ) 
NXX(I) 


CARO  3 


11-15  TIME  axis  scale  CODE  *iTm  FOLLOWING  SIGNIFI¬ 
CANCE  — 

isCxd)  so  automatic  scale  for  time  axis. 

«1  SCALE  TIME  AXIS  ACCORDING  TO  SPEC- 
IFiFn  PARAMETERS. 

16-25  minimum  VALUE  Of  TIMF. 

26-35  maximum  value  or  TIME. 

36-40  SCALE  FACTOR  (In  TFRmS  of  10**NXX>  TO  BE  USED  FOR 
THE  TIME  AXIS. 

(SEE  NOTE  (5)  FOR  CAGE  MOTION  PLOTS) 


ISCYU.I)  11-15 


riYd.n  if>-25 
YFY(l.I)  26-35 
NYY(l.I)  36-AO 


SCALE  CODE  FOR  THE  Y  AXIS  NO.  1  WITH  THE  FOLLOW¬ 
ING  SIGNIFICANCE  — 

I SCy  !  1  *  I > «0  AUTOMATIC  SCALE  FOR  Y  AXIS  NO.  1. 

■1  SCALE  Y-AXIS  NO. I  AS  SPECIFIED  And 
SKIP  AXIS  IF  data  OUT  OF  SCALE. 

■2  SCALE  Y-AXIS  N0.1  AS  SPECIFIED  AmD 
PLOT  FVEN  IE  DATA  OUT  OF  SCALE. 

■3  SCALE  Y-AXIS  NO, 1  AS  SPECIFIED  and 
TRUNCATE  DATA  IE  OUT  OF  SCALE. 
MINIMUM  VALUE  OF  the  VARIABLE  TO  RE  PLOTTED  On  Y 
AXIS  NO.  1. 


maximum  VALUE  Of  TnF  VARIABLE  TO  BE  PLOTTED  On  Y 
AXIS  NO.  1. 

SCALE  FACTOR  <!'•  TfR'-S  OF  10**NXx)  TO  RE  USED  FOR 
THE  X  AXIS. 


CARDS  4  AND  5 


RDREBP59 

RDREBP60 

RDRERP61 

RDREBP62 

RDREBP63 

R0REBP64 

PDRERP65 

RDREBP66 

RDRERPF7 

RDRERP68 

RDRERP69 

RDRERP70 

RDREBP71 

RDRERP72 

R0»EBP73 

RDREPP74 

RDRERP75 

RDRERP76 

R0RERP77 

RDREBP78 

RDRERP79 

RDRERP80 

RDRERPA1 

RDREBP82 

RDRERPR3 

RDBEPP°4 

RDRERPP5 

RD»EnpP6 

PDBErppt 

RDREBPP8 

RDRERPP9 

rdrerpro 

RDPERP9] 
RDREBP92 
RDRERP93 
RDPEBP94 
RDRERPR5 
RDRERP96 
RDREBp97 
RDRERP98 
RDREPP99 
RDREB1 DO 
RDRER101 
ROREB1 02 
RDREB103 
ROPEB104 
RDRER10S 
RDREB106 
RORER107 
RDREB108 
RDREB109 
ROPER  1 1 0 
RDREall  1 
RDpEpn2 
RD»EB1 13 
RDRER 1 1 4 
RDREB115 
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C  — — - - - -  RDREB 1 16 

C  SIMILAR  TO  CAPO  3  TOP  v  AXIS  NUMBERS  2  ANO  3.  RDRER1 1 ? 
C  THE  Y  AXES  ARE  NUMBERED  FROM  THE  ROTTOm  TO  THE  TOP  OF  THE  PLOT,  R0PEB118 
C  '  RDREB1 1 9 
C  ROREB120 
C  DATA  FILE  MANAGEMENT  RDREB121 
C  »Mii>i»«tnnni>  RDREB122 
C  ROREB123 
C  THE  FOLLOWING  DATA  FilES  ARE  USED  DURING  EXECUTION  OF  DREBP  —  RDREB 124 
C  RDREB125 
C  SOL  -  the  MASTER  Data  FILE  CONTAINING  ALL  THE  SOLUTIONS  for  THE  RDREB126 
C  GIVEN  BEARING  ELEMENT.  RDREB127 
C  SOLI  RDRE8128 
C  TO  -  LOCAL  FILES  CONTAINING  RESPECTIVE  INFORMATION  FOR  PLOT  RDREB129 
C  SOL  10  NO.  1  TO  10  RDREB130 
C  RDREB131 
C  A  MAXIMUM  OF  TEN  PLOTS  IS  PRODUCEO  BY  -ROPEBP-,  THE  ACTUAL  NUMBER  ROREB132 
C  OF  PLOTS  MAY.  HOWFVER.  DEPEND  ON  THE  PARTICULAR  MASTER  RDREB133 
C  DATA  FILE  SOL  WHICH  IDENTIFIES  THE  BEARING  ELEMENT  AND  HENCE  ThE  PDREP1 34 
C  REQIREO  PLOTS.  '  RDREB135 
C  RDREB136 
C  OUTPUT  DESCRIPTION  RDREB137 
C  RDREB138 
C  RDREBl 39 
C  THE  PRIMARY  OUTPUT  Or  ORERP  CONSISTS  OF  THE  VARIOUS  PLOTS  FOR  THE  RDREB140 
C  DIFFERENT  BEARING  PARAMETERS.  AS  DISCUSSED  ABOVE  THE  BEARING  RDREB1M 
C  PARAMETERS  ANO  ThF  DUMBER  OF  plots  ARE  DETERMINED  from  The  DAIa  RDRE8142 
C  StlPPLlFO  BY  The  MASTER  DATA  FILE  SOL.  IN  GENERAL  FOR  THE  VARIOUS  RDREB 1*3 
C  DATA  FILES  GEnEPaTEU  BY  DpEB*  THE  number  of  PLOTS  GENERATED  BY  RDRER144 
C  DREBP  ARE  INDICATED  BELOW  —  RDRER145 
C  RDREP1A6 
C  MOTION  OE  ROLLING  ELEMENT  10  PLOTS  RDREB 1 47 
C  MOTION  OF  CAGE  9  PLOTS  RDREB 1 48 
C  MOTION  OF  OUTER  RACt  A  PLOTS  RDREB 1 49 
C  MOTION  OF  INNER  RACE  6  PLOTS  RDREB150 
C  OVERALL  BEARING  PARAMETERS  1  PLOT  RDRE8151 
C  RDREB 152 
C  EACH  SET  OF  PLOTS  H*S  A  TITLE  PAGE  WHICH  CONTAINS  ThE  BEARING  IN-  RDREB 153 
C  FORMATION  AND  ALL  T«E  SCALE  FACTORS  FOR  THE  VARIOUS  PARAMETERS,  RDREB154 
C  RDREB 155 

c  along  with  all  thf  Rlots  orebp  also  produces  a  brief  print  output.  rdrebi56 
c  a  listing  of  all  data  cards  as  supplied  by  the  user  is  printed,  rdrebi57 

C  FOR  DIAGNOSTIC  PURPOSES.  AS  THE  CARDS  ARE  READ  IN.  THE  ESSENTIAL  RDREB158 
C  PRINT  OUTPUT  CONSISTS  OF  THE  BEARING  IDENTIFICATION.  FUNDAMENTAL  RDREB159 
C  SCALE  FACTORS  ANO  A  OETAlLED  TABLE  CONTAINING  ALL  ThE  SCALE  IN-  RDREP160 
C  FORMATION  FOR  THE  VARIOUS  PLOTS.  THE  EXECUTION  OF  A  PRINT  STATE-  RDREB161 
C  MENT  CORRESPONDING  To  A  GIVEN  PLOT  NUMBER  ALSO  INDICATES  A  SUC-  PDRER162 
c  cessful  plot.  rdrebi63 

c  RDREB 164 
C  EXTERNALS  CALLED  RDREB 165 
C  »fi(t»tiin«n  RDREB166 
C  RDREB167 
C  OTHER  than  The  CONVENTIONAL  CALCOMP  SOFTWARE.  THE  FOLLOWING  ex-  RDREB 168 
C  TEPNAL  ROUTINFS  muST  PE  SUPPLIED.  T*F  LIST  INCLUDES  ThE  ROUTINES  RDREB169 
C  CALLED  DIRECTLY  By  URpBP  OR  Th£  ONE*  CALLFO  BY  ANY  ROUTINE  INIT-  RDREP170 
C  IALLY  CALLED  BY  OREdP.  RDRE°171 
C  ROREB172 


-108- 


ooooooooooooonoooooooooonoononoooooo  on oo  no  oooooononoo n  o  o  o  o 


c 


1.  AXC 

2.  FNUM0R 

3.  MINMAX 

4.  MTI 

5.  SFAC 

6.  SIGF 

7.  TITLE 


COMPUTE  SCALE  FACTORS  FOR  GIVEN  AXIS. 

PLOT  REAL  NUMBER  IN  CONVENTIONAL  E  FORMAT. 

find  minimum  and  maximum  values  for  a  given  axis. 

PLOT  MTI  INSIGNIA. 

PLOT  SCALE  FACTOR  ON  A  GIVEN  AXIS. 

FXtRACT  SIGNIFICANT  FIGURES  OF  A  REAL  NUMBER. 
SET-UP  ALL  PLOT  TITLFS. 


notes 


<1)  A  MAXIMUM  OF  3000  POINTS  may  pr  PLOTTED  IN  A  GIVEN  CURVE. 

IF  THE  NUMRFR  OF  POINTS  .GT.  3000  THEN  EITHER  THE  DIMENSION 
STATEMENTS  IN  THE  UNLABELLEU  COMMON  (APPEARING  IN  THE  MAIN 
PROGRAM  and  Trt£  SUBROUTINE  -«Tnmax->  MUST  BE  CHANGED  OR  THE 
START.  FINISH  AND  NUMBER  OF  POINTS  TO  BE  SKIPPED  SHOULD  BE 
ADJUSTED  (SEE  CARD  1)  SO  THAT  NO  mope  THAN  3000  POINTS  ARE 
PLOTTEO. 

(2)  APPROPRIATE  VALUES  FOR  -MACH-,  -TAB-  AND  -SMAR-  MUST  BE 
SET  IN  THE  OATa  STATEMENT  FOR  PPOPEP  PROGRAM  OPERATION  — 

MACH  3  1,  Ten  characters  per  word  for  the  CDC  SYSTEM. 
mach  »  2,  four  characters  per  word  for  the  ibm  system. 

TAB  a  WIDTH  (INCH)  OF  PAPER  BETWEEN  TEAR  LINES  U 
ANY.  TAB  MUST  °E  .  C-E.8.5.  ALL  PLOTS  ARE 

designed  to  fit  ok  *  8.5  inch  wide  paper. 

SMAR  S  •  width  of  margin  iisrn  BY  THE  SYSTEM  FOR  USE* 

Identification*  fop  no  system  margin.  smaw=o. 

(3)  THE  FIRST  FOUR  COMPONENTS  OF  ARRAY  -IFOR-,  DEFINED  IN  A 
DATA  STATEMFNT,  APE  USEO  0»"LY  WITH  MACH«1  FOR  THE  CDC 
SYSTEM.  IF  THE  10H  FIELD  IN  ThE  DATA  STATEMENT  IS  UNACCEPT¬ 
ABLE  ON  the  IdM  SYSTEM  WITH  maCh*2  THEN  THE  VALUES  FOR  The 
FIRST  FOUP  COMPONENTS  MAY  RF  trumcaTED  BUT  UNDER  NO  CIRCUM¬ 
STANCE  the  dimensions  of  -IFOR-  may  be  altered. 

(4)  THE  FORM  OF  THE  CALCOMP  INITI ALlZATIONN  ROUTINE  -PLOTS-  MUST 
BE  CHECKED  FOR  the  GIVEN  INSTALLATION  AND  THE  STATEMENT  72 
IN  -ROREBP-  MUST  BE  AOJUSTED  ACCORDINGLY. 

(5)  PLOT  NO.  9  FOR  THE  CAGE  MOTION  IS  A  CAGE  MASS  CENTER  ORBIT 
PLOT  AND  IT  CONTAINS  ONLY  TWO  AXES  CORRESPONDING  TO  THE  Y 
ANO  Z  POSITION  OF  THE  CAGE  MASS  CENTER.  IF  ANY  SCALING 
INFORMATION  is  SUPPLIED  FOR  This  PLOT  THEN  THE  DATA  FO«  Y 
POSITION  SHOULD  BE  SUPPLIED  Om  CARO  2  AND  THE  DATA  FOR  l 
position  should  be  supplied  om  capo  3*  which  is  normally 
USED  FOP  Y  axis  NO.  1.  ALSO  Jn  THIS  CASE  ISCY(l.I)  MAY  BE 
EQUAL  TO  ONLY  A  0  OR  1,  ThE  OTHFP  TWO  OPTIONS  ARE  NOT 
AVAILABLE  FOR  THIS  PLOT. 

(A)  -ROREBP-  REOUlRES  A  CENTRAL  mfmOoy  OF  ABOUT  164000B  WORDS 
ON  A  CDC  SYSTEM. 
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APPENDIX  B 


TYPICAL  COMPUTER  OUTPUT  FOR  THE 


100  mm  AND  150  mm  DMA  BEARINGS 
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NET  BR6  LOSS  <N*M/S>  *  b.lB?F*03  NET  L0AD*Sl!P  <N»M/S»  ■  9.324E«04 
CURRENT  LIFE  (HOURS)  =  1.536F»0J  OUTER  RACE  FIT  <H)  »  l.OOOE-OS 
INTERNAL  CLEARANCE  (M>  *  1.067E-04  INNER  RACE  FIT  (M)  *  9.7b7E-07 
O.R.  HOOP  <N/M»»2)  *  -1.427E*07  I.R.  HOOP  (N/M*»2)  »  l,l«3E»08 
CAGE  HOOP  (N/M**2)  *  3,SllF»07 


